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We present an algebraic framework for interacting extended quantum systems to study complex phenomena 
characterized by the coexistence and competition of different states of matter We start by showing how to 
connect different (spin-particle-gauge) languages by means of exact mappings (isomorphisms) that we name 
dictionaries and prove a fundamental theorem establishing when two arbitrary languages can be connected. 
These mappings serve to unravel symmetries which are hidden in one representation but become manifest in 
another In addition, we establish a formal link between seemingly unrelated physical phenomena by changing 
the language of our model description. This link leads to the idea of universality or equivalence. Moreover, 
we introduce the novel concept of emergent symmetry as another symmetry guiding principle. By introducing 
the notion of hierarchical languages, we determine the quantum phase diagram of lattice models (previously 
unsolved) and unveil hidden order parameters to explore new states of matter Hierarchical languages also 
constitute an essential tool to provide a unified description of phases which compete and coexist. Overall, our 
framework provides a simple and systematic methodology to predict and discover new kinds of orders. Another 
aspect exploited by the present formalism is the relation between condensed matter and lattice gauge theories 
through quantum link models. We conclude discussing applications of these dictionaries to the area of quantum 
information and computation with emphasis in building new models of computation and quantum programming 
languages. 
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I. INTRODUCTION 

Unveiling the fundamental principles behind complex be- 
havior in matter is a problem at the frontiers of condensed 
matter physics and embraces cases like the high-temperature 
superconductors, heavy fermions, and low-dimensional elec- 
tron liquids. Generically, describing the structure and behav- 
ior of matter involves studying systems of interacting quantum 
constituents (bosons, fermions, spins, gauge fields) whose 
fields satisfy the basic laws of quantum mechanics. Never- 
theless, the plethora of complex phenomena exhibited by na- 
ture exceeds our ability to explain them, in part, because the 
whole is not necessarily the sum of its parts [l|] and thus typi- 
cal perturbation-like-theory arguments or standard mathemat- 
ical techniques are not appropriate to disentangle its myster- 
ies. In this paper we present a unifying algebraic framework 
for interacting extended quantum systems that enables one to 
study complex phenomena characterized by the coexistence 
and competition of various states of matter. 

The emergence of such complex phenomena may be the 
result of very simple, undiscovered, principles that conspire 
against any straightforward explanation. Nonetheless, we ex- 
pect that two pillars of modern science, symmetry and topol- 
ogy, are key guiding principles behind those fundamental laws 
of emergence. Group theory and geometry have been funda- 
mental to the physics of the twentieth century and we count 
on them to continue playing such a role. Indeed, the notion 
of symmetry and its breakings have shaped our current con- 
ception of nature. Understanding the idea of invariance and 
its corresponding conservation laws turns out to be as funda- 
mental as determining the causes that prevent such harmony 
and leads to more complex behavior. Other kinds of order, 
not described by broken symmetries, are possible. For ex- 
ample, quantum orders of topological nature distinguish the 



internal structure of the state without breaking any local sym- 
metry, and its study is beyond the scope of the present paper 
0]. Another example is provided by the notion of emergent 
symmetries where, as the name indicates, the number of sym- 
metries of the state increases as the temperature is lowered, 
contrary to the broken symmetry case. These concepts are at 
the heart of the field of quantum phase transitions that studies 
the changes that may occur in the macroscopic properties of 
matter at zero temperature (i.e., T = 0) due to changes in the 
parameters characterizing the system. 

In this regard, the development of exact algebraic methods 
is one of the most elegant and promising tools toward the com- 
plete understanding of quantum phases of matter and their cor- 
responding phase transitions. Typically, there are no distinct 
length or time scales which separates the different competing 
orders and often these systems are near quantum criticality 
which makes their study extremely complicated, if not impos- 
sible, by the traditional techniques, such as standard mean- 
field or perturbation theories. The reason which prevents the 
effective use of these theories is precisely the key for the suc- 
cessful application of the algebraic methods, namely, the ab- 
sence of a small parameter for the various complex quantum 
orderings. In other words, one cannot systematically apply 
renormalization group ideas and easily integrate out irrele- 
vant degrees of freedom. On the other hand, those competing 
orders are frequently related by symmetry principles charac- 
terizing the critical behavior caused by the competing inter- 
actions, thus increasing the symmetry group of the effective 
low-energy physics. 

In the present manuscript we will be concerned with quan- 
tum lattice systems. A quantum lattice is a composite system 
identified with where Ng is the total number of lattice 
sites, and associated to each lattice site (or mode) i e iJ^' is 
a Hilbert space Jii of finite dimension D describing the local 
modes. The total Hilbert space, i.e., state space of the lat- 
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tice, is the tensor product of the local modes state spaces, 
Ti, = (g). "Hi, in the case of distinguishable subsystems. (No- 
tice that TC may support inequivalent tensor product decom- 
positions.) Indistinguishability places an additional constraint 
on the space of admissible states which manifests itself in their 
symmetry properties (e.g., antisymmetric states in the case of 
fermions) and, consequently, the physical Hilbert space H is 
a subspace of Hi- A pure state of the system (at T = 0) 
is simply a vector | ^E*) in H, and an observable is a selfadjoint 
operator O : H H. The dynamical evolution of the system 
is determined by its Hamiltonian H. The topology of the lat- 
tice, determined by the connectivity of the graph induced by 
the interactions in H, is an important element in establishing 
complexity. Unless specified, we will always consider regular 
lattices of coordination z. In the case of quantum continuous 
systems we can still use the present formalism after discretiz- 
ing, for example, the space coordinates. Going beyond this 
approach is deferred for a possible later publication. 

Every given physical system is represented by a language 
which consists of a set of operators that acts irreducibly on 
the Hilbert space H. For example, if our system consists of 
a collection of interacting electrons a natural language is a 
set of spin 1/2 fermion field operators. Defining a language 
amounts to establishing the state space H used to describe 
the physical system, while a Hamiltonian written in terms of 
that language comes to setting up the quantum dynamics ex- 
pected to depict the phenomenon of interest. Can we connect 
the different (spin-particle-gauge) languages of nature within 
a single algebraic framework? The answer to this question is 
yes and the key outcome is a set of dictionaries relating the 
languages representing the different physical systems. More 
precisely, by dictionary we mean an isomorphic (one-to-one) 
mapping connecting two languages. The existence of dic- 
tionaries provides not only a tool to explore complexity but 
leads naturally to the fundamental concept of universality, 
or equivalence, meaning that different physical systems dis- 
play the same behavior. The concepts of language and dic- 
tionary are introduced in section |ll] Previous to our work 
there were two seemingly unrelated examples of these types 
of mappings: The Jordan- Wigner (1928) U and Matsubara- 
Matsuda (1956) ||4|| transformations. We have not only gener- 
alized these (SU (2)) transformations to any spin, spatial di- 
mension and particle statistics but have proved a fundamen- 
tal theorem that permits the connection of the generators of 
the different languages which can be used to describe a given 
physical problem. 

SectionlUstarts defining the concept of a bosonic language. 
Then we prove our fundamental theorem and establish the 
necessary conditions for application of a given bosonic lan- 
guage to describe the physical system under consideration. 
The main result of this section is a proof of the existence 
of mappings between the generators of any pair of languages 
which belong to the same class, i.e., that can be used to de- 
scribe the same physical system. The necessary and sufficient 
condition for any pair of languages to be connected is that the 
dimension of their local Hilbert spaces D be the same. In ad- 
dition we prove that for each class of languages, there is at 
least one that is realized by the generators of a Lie algebra. To 



complete this picture, we introduce the notion of hierarchical 
language and prove that there is at least one hierarchical lan- 
guage in each class. Those languages are, in general, the most 
convenient ones to characterize the different quantum phases 
of the system under consideration. 

In section Uni we extend these conclusions to non-bosonic 
systems by introducing operators which transmute the modes 
statistics. These operators have a local and a non-local compo- 
nent. The non-local part is a trivial generalization of the trans- 
formation introduced by Jordan and Wigner fsl] to map spins 
5=1/2 into spinless fermions. The local component trans- 
mutes the statistics associated to the interchange of particles 
which are sharing the same lattice site i. Adding the trans- 
mutation of statistics to the fundamental theorem completes 
the characterization of each class of equivalent languages. In 
other words, the classes of bosonic languages defined in sec- 
tion |ll| are expanded in section |llll to include fermionic, any- 
onic, or hybrid (para-) languages. 

The most natural choice of local Hilbert space Tii is the 
one spanned by a single site (or mode) i basis of dimension 
D. The possibility to decompose the total Hilbert space of the 
problem 7i into subspaces or subsystems which are not nec- 
essarily single sites opens up the possibility to generate other 
dictionaries with unforeseen applications. (Remember that Ti 
may support different tensor product decompositions.) Sec- 
tion |W] expands on this concept and show the simple case of 
a dictionary that uses as local Hilbert space a bond (i, j) state 
space Hi^T-ij = Hij of dimension D'^ {H = Wij), 
and which is mapped onto a site Hilbert space of the same 
dimension. 

SectionlVlis devoted to show explicit connections between 
equivalent languages. The main purpose of this section is to 
illustrate through examples the application of the dictionar- 
ies developed in the previous two sections. Out of the many 
possible transformations predicted by our fundamental theo- 
rem and the transmutation of statistics, we selected a few of 
them which are useful for the applications described in the 
later sections. However, it is important to note that the proce- 
dure described to find the mapping connecting two languages 
in the same class is very straightforward. In addition, we show 
that the fractional exclusion statistics algebras emerge natu- 
rally from the present formalism. 

Given a model Hamiltonian operator representing a quan- 
tum system, we can use our dictionaries to translate it into an- 
other equivalent Hamiltonian written in a different language. 
In other words, we can write down the same Hamiltonian op- 
erator in many distinct ways according to the different lan- 
guages that are included in the corresponding class. Since 
some of these languages are naturally associated to specific 
physical entities (like spins, particles, gauge fields, etc.), the 
corresponding translation provides a rigorous connection be- 
tween seemingly unrelated physical systems. This is the main 
subject of section IVII where we show examples of com- 
pletely different quantum lattice systems described by the 
same model. Another interesting aspect of these mappings 
is the potential to unveil hidden symmetries of the Hamilto- 
nian, possibly leading to exact or quasi-exact solutions (sj. 
Sometimes particular languages allow us to recognize invari- 
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ant subspaces of our Hamiltonian. When the action of our 
Hamiltonian is restricted to an invariant subspace, the corre- 
sponding operator can always be written in a language which 
is more elementary than the original one (i.e., it realizes a re- 
duced number of degrees of freedom). In some cases, this pro- 
cedure enables one to recognize the hidden symmetries which 
lead to the quasi-exact solution of the model considered (i.e., 
the model is integrable within the invariant subspace). The 
quasi-exact solution of the one-dimensional t-Jz 0] model 
belongs to this class of problems. It is important to remark 
that this quasi-exact solution leads to the exact determination 
of the quantum phase diagram and the charge excitations of 
the t-Jz model. 

A new notion beyond Landau's concept of broken symme- 
try |0 01 is the complementary idea of emergent symmetry, 
i.e., the fact that new symmetries not realized in the Hamilto- 
nian describing the system can emerge at low energies. There 
are many instances in which the high-symmetry, low-energy, 
effective theory can be derived in an exact way, i.e., with- 
out appealing to approximate schemes like the renormaliza- 
tion group. For instance, the exact effective theory (ground 
state plus any charge excitations) of the one-dimensional t- 
Jz model is the S=\I2 XXZ model 0]. It is well known 
that the XXZ model has an infinite number of symmetries 
which make it exactly solvable by the Bethe ansatz. Another 
example is provided by the family of spin Hamiltonians for 
which the ground state is a product of spin singlets. This fam- 
ily includes the Majumdar-Ghosh model and many other 
generalizations These and other simple 

examples introduced in section lvTll illustrate the fundamental 
concept of emergent symmetry which, as we will show, can 
be used as a guiding principle to find new states of matter 
In a sense that will become clear later, this concept not only 
formalizes but also provides a systematics to the principle of 
adiabatic continuity that P.W. Anderson (11] has advocated as 
one of the two most important principles of condensed mat- 
ter physics (together with the concept of broken symmetry). 
States like the Fermi liquid or the band insulator which do 
not correspond to any broken symmetry can be characterized 
by their corresponding emergent symmetries. In this way one 
can establish a formal connection between different phases. 
For instance, the emergent symmetry of the band insulator is 
analog to the emergent symmetry of the singlet dimer mag- 
netic states. In most cases, the notion of emergent symmetry 
is approximate; however, still provides a guiding principle to 
identify the relevant degrees of freedom, and the nature of the 
ground state and low-energy excitations. 

Coexistence and competition of different quantum order- 
ings associated to broken continuous symmetries is one of the 
main subjects of section IvTIll There we show the fundamen- 
tal role played by the hierarchical languages for the classifi- 
cation of the possible order parameters and the calculation of 
the quantum phase diagram of a given model. We note that the 
hierarchical languages are the most natural ones to provide a 
unified description of the order parameters characterizing each 
phase. To illustrate the procedure we consider the example of 
the Heisenberg SU (N) Hamiltonians written in different lan- 
guages. We then take the local SU (N) order parameter and 



reduce its components according to the different subgroups of 
SU{N) that can be used to generate an equivalent language 
(for instance SU{2)). As an illustration, we show that the lo- 
cal order parameter for an S=l SU (2)-spin Hamiltonian can 
be either the usual magnetization or a spin-nematic order pa- 
rameter To conclude this section, we show how to obtain 
the quantum phase diagram of the bilinear-biquadratic 5=1 
Heisenberg model (for spatial dimensions d > 1) just by writ- 
ing the Hamiltonian in the SU{3) hierarchical language. It 
is important to remark that this zero temperature phase dia- 
gram was only known for semi-classical spins fl4ll . The fact 
that there is no calculation involved in the elaboration of this 
quantum phase diagram permits the reader to appreciate the 
power and potential applications of this algebraic framework. 

The formal connection between lattice models in con- 
densed matter physics and lattice gauge theories of high- 
energy physics is described in section IIXI There we take 
advantage of the existing quantum link models and the con- 
nection between spins and other degrees of freedom, such as 
gauge fields, which emerges from our algebraic approach. 

In addition to the fascinating field of quantum phase tran- 
sitions and statistical mechanics, our algebraic approach can 
be applied to the field of quantum information ITsll . This 
new paradigm of information processing takes advantage of 
the fundamental laws of quantum mechanics to perform op- 
erations which can be done at least as efficiently (i.e., with 
polynomial complexity) as with classical devices. The device 
that performs the manipulation of information is named quan- 
tum computer and the basic unit of information is the qubit 
(i.e., a two-level system). The quantum computer (a quantum 
many-body system in disguise) consists of an isolated set of 
quantum degrees of freedom which can be controlled and ma- 
nipulated at the quantum level to perform the required opera- 
tions. The different languages associated with the description 
of these quantum degrees of freedom are possible program- 
ming languages and indeed may realize different models of 
quantum computation. The development of dictionaries con- 
necting these languages is relevant not only to improve the 
efficiency of a quantum computer, but also to model and sim- 
ulate different physical phenomena. The applications of our 
algebraic framework to the fields of quantum information and 
quantum computing are discussed in section ^ There we 
show how each physical realization of a quantum computer 
has a class of languages associated to the realization of the 
quantum operations. For instance, if we consider the typical 
case of interacting S'=l/2 spins the most natural language is 
the one generated by the Pauli matrices. Hence we can use 
the dictionary provided by the Jordan- Wigner (JW) 0] trans- 
formation to simulate a fermionic system ||, 1 ^ . In this regard, 
the purpose of this paper is to generalize this idea to other pos- 
sible realizations of quantum computers. We also discuss the 
potential applications of this algebraic framework to the de- 
scription of recent experiments on bosonic systems in optical 
lattices. 

Finally, in section PTTl we summarize the main concepts in- 
troduced in this paper 
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II. GENERAL LIE ALGEBRAIC FRAMEWORK 

The notions of Hilbert space and linear maps or operators 
are central to the formulation of quantum mechanics flTil . 
Pure states (wave functions) are elements of a Hilbert space 
and physical observables are associated to Hermitian (self- 
adjoint) operators which act on that space, and whose eigen- 
values are selected by a measuring apparatus. The role of lin- 
ear operators in quantum mechanics is not restricted to the 
representation of physical observables. Non-Hermitian oper- 
ators are very often used. That is the case, for instance, of 
the creation and annihilation operators in second quantization. 
Moreover, the unitary operator that describes the time evolu- 
tion of the system is in general non-Hermitian. These observa- 
tions simply tell us that linear operators in general (Hermitian 
and non-Hermitian) play a more general role in quantum me- 
chanics since they provide the mathematical language which 
is required to describe the dynamical behavior of our quantum 
system. 

What do we mean by the term mathematical language! In 
order to provide a rigorous answer for this question it is im- 
portant to first determine what is the mathematical structure 
associated to the set of quantum operators. Quantum opera- 
tors form a complex linear (vector) space under the sum and 
the multiplication by a scalar over the field of complex num- 
bers. If we augment this vector space with a bilinear opera- 
tion (product □ between two operators), the set form an al- 
gebra. Quantum mechanics also requires this operation being 
non-commutative and associative. This additional structure 
makes the set of quantum operators form an associative alge- 
bra. In principle, any operator in this algebra can play a role 
in the description of our quantum system. However, one can 
always select a subset (which forms a basis) in such a way 
that any quantum operator can be expressed as a function of 
the selected ones. For instance, we can use the Pauli matrices 

to express any quantum operator associated with a spin 1/2. 
The particular subset of operators that we choose to express 
any quantum operator is the mathematical language that we 
will use for the quantum description of our system. The ele- 
ments of this subset will be called generators of the language. 

What are the conditions a given set of quantum operators 
must satisfy to become a language? How many different lan- 
guages can be used to describe a quantum system? What is the 
connection between the different languages? What is the most 
appropriate language to describe a particular system? What 
is the relation between language and symmetry generators? 
How can languages help us understand the phenomenon of 
universality? (It turns out that the notions of language and 
universality are closely related whenever a common language 
can be used to describe seemingly unrelated physical phenom- 
ena.) A great part of this paper is dedicated to answer these 
questions. In particular, in this section we introduce the no- 
tion of bosonic language, we establish a formal connection or 
dictionary between the different bosonic languages associated 



with a given bosonic system, and we also establish a formal 
relation between bosonic languages and Lie algebras. 

One of the fundamental steps toward a deeper understand- 
ing of classical systems was the recognition that dynamical 
variables which define the phase space are generators of con- 
tinuous transformations. The set of continuous transforma- 
tions forms a group and the infinitesimal generators provide 
a basis for a Lie algebra which is related to the group. A 
group is a non-empty set which is closed under an associative 
product □, it contains an identity, 1, and all of its elements 
are invertible. When the transformations in the group leave 
the equations of motion invariant the set becomes a symmetry 
group for the considered system. A real (complex) Lie algebra 
£ is a linear space over the field F of real (complex) numbers 
which is closed under a non-associative Lie product [ , ] that 
satisfies (a, b, c G £ and a, /? G F): 

[aa + /?b,c] = a[a, c] +/3[b,c] 
[a, b] = -[b,a] 

= [a,[b,c]] + [b,[c,a]] + [c,[a,b]] (2) 

Whenever a continuous transformation is a symmetry of our 
physical system, the infinitesimal generator, which is an ele- 
ment of the corresponding Lie algebra, becomes a conserved 
quantity. As a consequence, it is convenient to choose that 
quantity as one of the coordinates for the phase space of our 
system. 

In the same way, it may be appropriate to choose the quan- 
tum operators (generators of the language) for the description 
of a quantum system in such a way that they simultaneously 
provide a basis for a Lie algebra. In other words, the use 
of symmetry generators for the description of an interacting 
quantum system incorporates symmetry as a guiding principle 
to find out the possible solutions. As we demonstrate below, 
it is always possible to find a language whose generators form 
a basis of a Lie algebra. 

We start by considering the bosonic languages of quantum 
mechanics. To introduce the definition of a bosonic language, 
we need to define first the concept of monoid. A monoid is 
a triple (9JI, □, 1) in which 971 is a non-empty set, □ is an 
associative product in Tl, and 1 is an element of 371 such that 
□ (1, a) = a = □(a, 1) for all a e Tl. In this way we see that 
the the concept of monoid generalizes the notion of g roup ; a 
group is a monoid all of whose elements are invertible lllal . A 
bosonic language is a set of operators which can be grouped in 
subsets Si (associated to each mode) and satisfy the following 
conditions: 

• Each element of Si (/i S [l,iVg]) belongs to the 
algebra of endomorphisms for the vector space Hi over 
the field of complex numbers C, : Hi ^ Hi, and 
these elements are linearly independent. 

• The elements of Si generate a monoid of linear trans- 
formations under the associative product in the algebra 
which acts irreducibly on Hi in the sense that the only 
subspaces stabilized by Si are Hi and (0 is the null 
vector). 
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• If and are elements of different subsets Si and Sj, 
then bfby = □(bf , bp = □(bj', bf ). 

The elements of the sets Si are called generators of the 
bosonic language. Combining the associative product and the 
additive operation in the algebra of endomorphisms of Hi, 
EndcTii, we can define the non-associative Lie product [ , J , 
which is called commutator: [bf,by] = b^bj - b^f (ij. 
Then, the last condition can be reformulated by establishing 
that the commutator of elements in different subsets is zero 

[b[',by] = 0, if i^j. (3) 

It is also important to notice that the set Si is not necessar- 
ily closed under the regular product (composition) or the Lie 
product (commutator). If the set Si is closed under the Lie 
product, the elements of Si generate a Lie algebra. We will 
denote the Lie algebra by Si. In addition, since each generator 
is represented by an endomorphism of Hi there is a particular 
representation of Si associated to the bosonic language. 
The second condition for a bosonic language implies that 
is irreducible. The third condition implies that the global Lie 
algebra associated to the complete set of generators is the di- 
rect sum of local algebras Si, S = 0; Si. Therefore, if the 
set Si is closed under the Lie product, we can represent the 
bosonic language by the conjunction of the Lie algebra S, and 
the irreducible representation : 5 A . The dimension of 
F5 is equal to the dimension of the local Hilbert space Hi'. 
dimF^ = D. The algebra of endomorphisms for the vector 
space Hi is the enveloping algebra of Si. 

The following fundamental theorem shows that two lan- 
guages are equivalent if they have in common the dimension 
D of their local Hilbert space Hi, and D is finite. In other 
words, there is always a dictionary connecting both languages, 
and a given physical phenomena can be described with either 
one. A corollary of the theorem is that given a particular lan- 
guage one can always find an equivalent one which can be 
expressed as 5 A F5, i.e., the operators are generators of a 
Lie algebra Si in the irreducible representation F5. In general 
there are various complex Lie algebras which can be used to 
build equivalent languages. The only condition is that they 
have to admit an irreducible representation of dimension D. 
Clearly, the number of Lie algebras satisfying this condition 
increases with D. 



A. Equivalent Classes of Bosonic Languages 

The demonstration of the fundamental theorem of this sec- 
tion is a direct consequence of the classical theorem of Burn- 
side I20I1 . which plays an important role in the theory of rings. 
Given the above definitions we can enunciate the theorem in 
the following way: 

Burnside 's theorem. Let G be a monoid of linear transfor- 
mations in a finite dimensional vector space V over an alge- 
braically closed field F, that acts irreducibly on V in the sense 
that the only subspaces stabilized by G are 1/ and 0. Then G 
contains a basis for EnApV over F. 



EnApV is an abbreviated notation for the ring of endomor- 
phisms of V over the field F. In quantum mechanical sys- 
tems, the field F corresponds to the complex numbers C. The 
demonstration of the Burnside's theorem can be found, for 
instance, in Ref. Ii20i1 . An immediate consequence of this the- 
orem is the following one which is the basis for connecting 
different bosonic languages: 

Fundamental theorem ("On the equivalence of languages"). 
Given two bosonic languages having the same finite dimen- 
sion D of their local Hilbert spaces. Hi, the generators of one 
of them can be written as a polynomial function of the gener- 
ators of the other language and vice versa. 

Proof: The proof is a trivial application of Burnside's the- 
orem. We need to notice first that since the dimension D is 
the same for both languages, the spaces Hi are also the same 
(vector space of dimension D over the field of complex num- 
bers (D). Let us consider the monoid of transformations Gi 
generated by multiplying the generators of the first language 
until the enlarged set becomes closed under the product. The 
second condition in the definition of a bosonic language states 
that Gi acts irreducibly on Hi. Since the dimension D of Hi 
is finite, Burnside's theorem guarantees that Gi admits a ba- 
sis for End(c7ii. Therefore, any endomorphism in Hi can be 
written as a linear combination of endomorphisms in Gi. In 
particular, the generators of the second language can be writ- 
ten in this way because they belong to End^Tii. Since each 
element of Gi is a product of generators of the first language, 
this concludes the demonstration. 

This theorem establishes an isomorphism between the al- 
gebras of endomorphisms associated to each of the two lan- 
guages. Motivated by this observation we can introduce the 
notion of classes of equivalent bosonic languages. We will 
say that two bosonic languages belong to the same class if 
they have the same dimension D of their local Hilbert spaces 
Hi. The fundamental theorem establishes the existence of dic- 
tionaries connecting languages within the same class. As a 
consequence, we can use any bosonic language in the class 
to describe a given physical phenomena. The natural ques- 
tion which emerges from this result is: What is the most ap- 
propriate language in a given class for describing our partic- 
ular problem? There is no generic answer to this question. 
Nonetheless, the following two corollaries give an important 
hint for problems which are invariant under particular trans- 
formations, because they relate the notion of language to the 
generators of symmetry groups. 

Corollary I: In each class of bosonic languages there is at 
least one which is the conjunction of a Lie algebra S and an 
irreducible representation F5 {S A F5), i.e., the generators of 
the bosonic language are generators of the Lie algebra Si in 
the representation F5. 

Proof. First, we need to notice that each class is character- 
ized by the dimension D of the local Hilbert space Hi. Let us 
consider the group U (1) ® SU (2). The Lie algebra associated 
to this group is Ci = u{\) su{2). The generators of Ci are 
{l;,Sf,Sy,S^i}, 

[S'f , Sj'] iSi^f-^j^uxSi , ^i,iy,X^ X, y, z (4) 

{S^ = Sf ± iS\ and e is the totally antisymmetric Levi- 
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Civita symbol), and there is one irreducible representation 
of dimension D for each possible value of D 111 lil . Therefore, 
the set Si = {1;, S'f , , S?} in the representation fulfills 
the three requirements for a bosonic language and operates in 
a local Hilbert space of dimension D. Since Si is the set of 
generators of the Lie algebra Ci, making iSi = Ci and = 
, we have proved the corollary I. 

The proof of corollary I shows that any bosonic problem 
with a local Hilbert space of dimension D can be described 
with S'[/(2)-spins of magnitude (representation) 5 = (£> — 
l)/2. The Matsubara-Matsuda (MM) t4'3 transformation is 
the simplest application of this corollary to D = 2. (Indeed, 
one can construct generalized MM transformations for any D 

m.) 

We introduce now another definition which is motivated by 
the next corollary. A given bosonic language will be called 
hierarchical if any local physical operator Oi can be written 
as a linear combination of the generators of the language, i.e., 

Oi^Y. ^M&r . (5) 

where € C, and it is the conjunction of a Lie algebra S 
and an irreducible representation Tg. 

Corollary II: In each class of bosonic languages there is one 
which is hierarchical and its generators are the identity and the 
generators of su{N = D) in the fundamental representation. 

Proof: For each class, with dimension D of the local Hilbert 
space Hi, we consider the group L/(1)(8)S'[/(A^) with = D. 
The generators of the Lie algebra associated to this group, 
Ci = u{l) Q su{N), are the identity plus the generators of 
su{N). Since the fundamental representation of Ci has 
dimension D ^ N, the conjunction of £ = 0j Ci and this 
representation is one possible language for the class consid- 
ered. Since the dimension of Ci is D^, which is the dimen- 
sion of the vector space End^Tii, then the generators of Ci 
also form a basis for EndcTii. 

The first consequence of corollary II is that the generators 
of any language can be expressed as a linear combination of 
generators of a hierarchical language in the same class. Again, 
the most trivial example is given by the class of bosonic lan- 
guages containing the spin 5=1/2 lattice. The generators of 
any language (like hard-core bosons or any two level system) 
in the same class can be written as a linear combination of the 
identity and the Pauli matrices. We will see later that corollary 
II is the key to get a hierarchical classification of the possible 
broken symmetries of a given physical Hamiltonian. 

We consider now two additional examples that illustrate in 
detail the contents of the fundamental theorem and the sub- 
sequent corollaries. The first example corresponds to hard- 
core bosons with Nf different flavors a. Since they are hard 
core only single occupancy is allowed, i.e., the eigenvalues of 
"i = "i" ^re either or 1 (blJ)}^^ = 0, and hia = blj>ia 
is the number operator for the flavor a at the site i). To distin- 
guish between site and flavor indices we will adopt the follow- 
ing convention: Sites are denoted by latin indices while the 
different flavors (or orbitals within the same site) are labeled 



by greek indices. The minimal set Si of operators that we can 
use to generate a bosonic language which is appropriate for 
hard-core bosons is; Si = {1;, bl^jb^^} with 1 < a < Nf. It 
can be shown that this set satisfies the three requirements for 
a bosonic language. The dimension of the local Hilbert space 
for these endomorphisms is D = Nf + 1. Then by varying 
the total number of flavors we can generate all possible values 
of D. Since each class of bosonic languages is characterized 
by the value of D, these hard-core bosons provide an exam- 
ple of a bosonic language in each class. It is clear that the 
set Si is not closed under the Lie product. Therefore, we can- 
not associate a Lie algebra to this minimal bosonic language. 
However, if we increase the number of generators in the set Si 
by including bilinear forms of the type fci^fc;^, then the new set 

= {hMaAcMJiph With l<a, f3<Nf, becomes 
closed under the Lie product 

' ^b^Ac^l = -niSap - bl^b^p) , (6) 

This means that the extended set Si is now a set of gener- 
ators for a Lie algebra in a particular representation. From 
the commutation relations (Eq. (|6j) we can conclude that Si 
is the direct sum of an algebra, generated by the iden- 
tity 1;, and an su{N) (N = D ^ Nf + 1) algebra gen- 
erated by {bl'bi^MJip}- Si = u{l)®su{N). The rep- 
resentation F^ is the fundamental representation of su{N) 
(dimF^ = N). Therefore, the new language is a hierarchi- 
cal one. In section IVAl we give a more detailed description 
of this particular language. Here we only want to emphasize 
the practical consequences of the fundamental theorem and 
its corollaries. The first non-trivial observation is that for each 
system of interacting multiflavored hard-core bosons, there is 
an equivalent system of interacting 5?7(A^)-spins in the fun- 
damental representation (with the minimal non-zero magni- 
tude). For Nf = 1 we recover the well-known MM trans- 
formation 0]. We can see now that this is the generalization 
to su{N) of the MM transformation. With this example we 
can envision the broad set of applications derived from the 
fundamental theorem. Another consequence of the theorem 
is that any physical theory for a bosonic system can be for- 
mulated in terms of multiflavored hard-core bosons if the di- 
mension of the local Hilbert space is finite. The usefulness 
of this formulation will depend on the particular system as it 
is illustrated in the next sections. Since the second language 
is hierarchical we can write down any local physical operator 
(endomorphism in Hi) as a linear combination of its gener- 
ators. Therefore, each of these generators will appear in the 
Hamiltonian under consideration with a power not larger than 
one. In addition, if the Hamiltonian has a global symmetry 
generated by a direct sum of local transformations: Ti, the 
symmetry will become explicit by writing the Hamiltonian in 
terms of a hierarchical language. The most basic example is 
the case of MM hard-core bosons (Nf = 1) in a lattice, de- 
scribed by a Hamiltonian with a kinetic energy term and a 
nearest-neighbor density-density interaction. The expression 
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for the Hamiltonian in terms of the first language defined by 
the set Si is 

iJxxz = tJ2iHh + blk) + VY.ini - - i) . (7) 

(iJ> (ij) 

where (i, j), refers to nearest-neighbors in an otherwise regu- 
lar d-dimensional lattice. Since 6; and 6; are not generators of 
a Lie algebra, the eventual global symmetries of ifxxz remain 
hidden in this particular language. However, if we translate 
Hxxz to the second 5J7(2)-spin language using the dictionary 
provided by Matsubara and Matsuda 0| 

Sr = 6; , (8) 

we can immediately unveil the hidden symmetries of iJxxz. 
The well-known expression for iJxxz in terms of the su{2) 
generators (i.e., the equivalent spin Hamiltonian) is 

H^ = J2 JzS?S^ + ^ {S^Sr + sr 5+) . (9) 

<iJ> 

The magnetic couplings, Jz and J±, are related to the original 
parameters, t and V, by the relations: = V and J± = 2t. 
It is clear from the last version of iJ>o<z (Eq. (|9}) that the orig- 
inal model has a global SU{2) invariance if V ~ 2t, i.e., 
it is in the isotropic Heisenberg point. The existence of this 
5?7(2)-symmetric point has a very important consequence for 
the phase diagram of the bosonic model of Eq. 0: If there 
is a charge density wave (CDW) instability at that point, the 
SU{2) invariance implies that there is also a Bose-Einstein 
condensation and vice versa. The order parameters of both 
phases are different components of a unique order parame- 
ter in the spin language, i.e., the staggered magnetization of 
the antiferromagnetic phase (t > 0). The z-component of 
the staggered magnetization is mapped onto the CDW order 
parameter for the bosonic gas, while the transverse compo- 
nent is equivalent to the order parameter for the Bose-Einstein 
condensation. Only one of these two phases, which are coex- 
isting at the SU{2) invariant point, is stable when we depart 
from the symmetric point in any of both directions in parame- 
ter space (Bose-Einstein condensation if ^ < 2< and CDW if 
V > 2t). In this very simple example we can see the advan- 
tages of using a hierarchical language {su{2) in this case). In 
the first place, we can immediately recognize the high sym- 
metry points. Secondly, we can describe an eventual broken 
symmetry state at those points in terms of a unified order pa- 
rameter 12311 . If we were to use a non-hierarchical language 
to describe the same problem, we would find coexistence of 
more than one phase at the high symmetry points. The order 
parameters of each of those phases are different components 
of the unified order parameter that we simply found with the 
hierarchical language. These ideas are developed in more de- 
tail in section IVlIII Here the aim is to give a flavor of the 
potential applications of our fundamental theorem. 



Using the concept of transmutation of statistics, to be in- 
troduced in the next section, we will extend the notion of 
classes of bosonic languages to more general classes contain- 
ing non-bosonic languages (i.e., languages for which Eq. Q 
is replaced by a different algebraic condition). The most well- 
known examples of non-bosonic languages are the fermionic 
ones for which the commutator of Eq. (|3} is replaced by an 
anticommutator We will see that there are non-local transfor- 
mations which map non-bosonic languages into bosonic ones. 
In this way, the simultaneous application of these transforma- 
tions and the fundamental theorem provides the natural path 
toward a unification of the different languages used to describe 
quantum systems. 



III. TRANSMUTATION OF STATISTICS 

In the previous section we have shown that two different 
bosonic languages can be connected (i.e., there is a dictio- 
nary connecting the two languages) if they have the same di- 
mension of their local Hilbert spaces D. However, we know 
that the bosonic languages do not exhaust all the possible lan- 
guages of quantum mechanics. The best known example is 
the fermionic language, whose creation and annihilation op- 
erators obey a closed set of anticommutation relations. There 
are many other examples of non-bosonic languages that we 
address below. 

Is it possible to connect these non-bosonic languages to the 
bosonic ones introduced in the previous section ? A positive 
answer to this question is given in this section by introduc- 
ing operators which transform commutators into anticommu- 
tators and vice versa. These operators have a local (related to 
the generalized Pauli exclusion principle (see section IV^ ) 
and a non-local (related to the exchange statistics) compo- 
nents. The local component is derived at the beginning of this 
section, while the non-local one is introduced in the second 
part. We will see that an Abelian anyonic statistics can also 
be achieved by introducing a continuous parameter in these 
transformations. (Non-Abelian statistics is beyond the scope 
of the present manuscript.) The transmutation of statistics to- 
gether with the fundamental theorem establish the framework 
needed to classify the classes of equivalent languages. This 
means that to describe a given physical problem one can use 
the original language or any other belonging to the same class. 

Before proceeding, we would like to mention that it seems 
as if there are two unrelated notions of anyonic statistics in 
the literature. One tied to quantum mechanics in first quan- 
tization in the coordinate representation, and another derived 
within the framework of quantum field theory. In both cases 
the original motivation to introduce such particles was basi- 
cally as an inherent possibility in the kinematics of (2h-1)- 
dimensional quantum mechanics and clearly the concepts, if 
correctly implemented, should be equivalent but are not in- 
deed I25II . In this manuscript, our anyon notion is consistent 
with the one developed in quantum field theory where the ex- 
clusion properties are preserved under statistical transmuta- 
tion. In this way, fermions can be kinematically transformed 
into hard-core bosons but not into canonical ones. 
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A. Fermionic Languages 

In addition to bosons, the other type of fundamental parti- 
cles found in Nature are the fermions. We have seen that the 
notion of bosonic languages is closely related to the concept of 
Lie algebras, since for each class of bosonic languages there 
is at least one language whose elements are generators of a 
Lie algebra. However, the same cannot be done for fermions. 
This can be easily seen by noticing that the third condition 
for a bosonic language (see Eq. Q) is not valid in general for 
fermions. In addition, the main consequence of the third con- 
dition for a bosonic language is that the global Lie algebra is 
the direct sum of the local ones associated to each subset Si. 
The generalization of these concepts to fermionic languages 
can be done by introducing the notion of Lie superalgebras 
(see for instance Ref. 12611 ). The fermionic languages are 
associated to Lie superalgebras in the same way the bosonic 
languages are associated to Lie algebras. Therefore, to give a 
general definition of a fermionic language we should first in- 
troduce and explain the notions of Grassman algebras, which 
are associative superalgebras, and Lie superalgebras. Since 
this is beyond the scope of the present paper, we will only 
consider the fermionic language generated by the canonical 
creation and annihilation operators 



(10) 



and other languages obtained by imposing particular local 
constraints on the canonical fermions. These generators, to- 
gether with the identity, generate the Heisenberg Lie super- 
algebra. In analogy to the bosonic languages (see Eq. Q), 
the Lie product (in the superalgebra) of two elements acting 
in different sites (different subsets Si,Sj ) is zero. Thus, in- 
stead of having a direct sum of local Lie algebras like in the 
bosonic case, we have a direct sum of local Lie superalgebras. 
In the case of canonical fermions the local Lie superalgebras 
are Heisenberg superalgebras. 



B. Local Transmutation 

We will start by considering multiflavor fermions C;^ (a £ 
[l,Nf]) which satisfy the canonical anticommutation rela- 
tions (Eq. ilQ\ ). Other types of fermions, usually considered 
in physics, can be derived from the canonical ones by impos- 
ing particular constraints. For this reason, the transformations 
derived for canonical fermions can be extended to these other 
fermionic algebras by incorporating those constraints. This 
procedure is illustrated with different examples in the next 
section. 

The canonical fermions can be transformed into bosons 
which are hard-core in each flavor (the eigenvalues of 

nia ~ ^L^iQ either or 1), i.e., two or more bosons of the 
same flavor are not allowed on the same site, but two of dif- 
ferent flavors are permitted. In the previous section we have 
shown that a physical theory for objects obeying commuta- 
tion relations (Lie brackets) can be formulated in terms of a 



bosonic language. By the present connection we will be able 
to extend this statement to fermions (or anyons, in general) 
through a transmutation of statistics. To this end, we need to 
define a different type of multiflavored hard-core boson 
satisfying the following commutation relations 

UhAp^ = [bLbU^o, ^^^^ 

which implies {bi^,bl^} = 1. The Lie algebra generated 
by these bosonic operators is £ = j su{2), i.e. each set 

{^L' bia^fiia — 1/2} generates an su{2) algebra. 

To show that these bosons (hard-core in each flavor) can be 
expressed as a function of canonical fermions, one defines 



(12) 



where "Tj^ = exp[i7r fijp] is the local transmutator, and 

/3<Q 

we are assuming a particular ordering for the flavor index a. 
From the expression for Tj^ it is clear that 



T- ^ I ^T- 



(13) 



It is straightforward to verify that the cj -operators satisfy local 
canonical anticommutation relations (Eq. ( I10> when i = j). 

In this way we have established a mapping between 
fermions and bosons which are operating locally (on a given 
orbital or mode j). In other words, we have related the subset 
'^j ~ {^ja'^ja'^ja ~ 1/2} of local generators of a bosonic 
language to the subset 5j = {cjo,, Cj^.^, nja — 1/2} of local 
generators of canonical fermions (njc = Cj^Cj^,). 



C. Non-Local Transmutation 

So far, we have only transmuted the commutation relations 
between generators which belong to the same site or subset 
S\. For commutation relations of two generators of different 
sites we need to introduce a non-local operator A'j. Jordan 
and Wigner Q] were the first to introduce such an operator 
in connection to their one-dimensional (Id or d = 1) trans- 
formation between spins S=\I2 and spinless fermions. The 
so-called kink operator that they introduced is 



A'j'* = exp[z7r ^ 

Kj 



"I 



(14) 



where n\ is the number operator for spinless fermions at the 
site 1. It is clear that for multiflavor canonical fermions we 
only need to replace n\\iy n\ — n\a, where a denotes the 
flavor Therefore, Eq. il2i must be replaced by 

ct = p. f.^ K^. = 6! /Ct , (15) 

where ifj = Kj'^ for a Id lattice. Even though this is a non- 
local operator, it does not introduce long-range interactions if 
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the model has only short-range terms |27]. In other words, 
a given Id Hamiltonian can be written in terms of bosons, 
fermions or anyons, and the interactions remain short range 
for all the cases. This is a special characteristic of one di- 
mension. The only consequence of changing the statistics of 
the particles is a change of the short-range interactions in the 
original basis. Therefore, the concept of particle statistics in 
one dimension becomes irrelevant since any physical system 
can be described with a bosonic language without changing 
the short-range character of the interactions lEsIl . 

The next step is the generalization of K;to higher dimen- 
sions. This has been done by Fradkin 12911 and Y.R. Wang 
llsoll . who considered the generalization of the traditional JW 
transformation for the two-dimensional {2d) case, and Huerta- 
Zanelli llsill and S. Wang ll32il . who did the same for higher 
dimensions. We will see that, in the same way we did for the 
Id case, these generalizations can be extended to transform 
canonical fermions into bosons and vice versa. The general- 
ization given by Fradkin L29.1 for the 2d case is 

A-f ==exp[i^a(lj)ni] , (16) 
1 

Here, a(l, j) is the angle between the spatial vector 1 — j and 
a fixed direction on the lattice, and a(j,j) is defined to be 
zero (see Fig. Again for the case of multiflavor canoni- 
cal fermions we just need to replace n\ by n\. We comment 
that the Id kink (or string) operators constitute a particular 
case of Eq. il6\ with a(l, j) = n when 1 < j and equals zero 
otherwise. 























1 










< 



















FIG. 1: Schematics of the geometry and function describing one of 
the possible ways of defining the statistical transmutator K?'' in 2d. 
The dotted line represents a fixed direction on the lattice. 

A general procedure for finding the string operators in d 
dimension has been given by Wang ll32il in the way we de- 
scribe below. Following the Id and 2d cases, the expression 
proposed for the operator Kj is 

Aj exp[i^w(lj) ni] , (17) 
1 
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where is a function to be determined by imposing the 

transmutation of statistics. It can be shown that this is equiva- 
lent to the antisymmetric condition 

exp[iw(lj)] = -cxp[iw(j,l)], if 1 ^ j 

w(l,l) = 0. (18) 

The simplest solution for a Id lattice is 

^(ij) = ^e(j-i) (19) 

where Q{x) is the Id step or Heaviside function. This is the 
solution found by Jordan and Wigner (see Eq. il4\ ). The so- 
lution introduced in Ref. ^9ll for the 2d case corresponds to 

t^(lj) =a(lj) . (20) 

However, it has been pointed out by Wang ll32il that this is not 
the only possible solution since 

= 7r[Qin - h){l - Si,,,) + - l2)5i,,,l (21) 

with 1 = liGi + ^262 and j = jiei + ^262, also satisfies 
Eqs. (I18> . The advantage of this solution is that its general- 
ization to higher dimensions is straightforward (1 = ^^e^, 
and jjL e [1, d\). For instance, in 3d we have 

+ Q{]z-h)KnKj.]- (22) 

In the next section we will see that, for d > 1, the operator 
ifj introduces non-local interactions in the new representa- 
tion of the original Hamiltonian. A phase factor of the form 
exp[i A • ds] appears in between the product cl^c^ of the 
kinetic energy-like terms. The field operator A(x) is defined 
by 

A(x) =^Va;(x,l)ni, (23) 
1 

where V represents the lattice gradient (/i, v G [1, d]) 

V^u;(x,l) =cj(x + e^,l) -tj(x,l). (24) 

It was shown ll29il that A(x) is the vector potential of a gen- 
eralized Chern-Simons construction for a lattice. The field 
strength associated to this vector potential is the lattice 
rotor 

F^,(x) = V^A,(x)- V,A^(x), 

Fo^(x) = VoA^(x)- V^Ao(x), (25) 

where Vo is the time derivative and Ao(x) is a scalar field. 
Since is a Lagrange multiplier field for the Chern-Simons 
Lagrangian, it can be integrated out to get the Gauss law 

fix = -Fi2(l) . (26) 

TT 
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This relation imposes the constraint between charge and flux 
giving rise to the vector potential A. The physical interpreta- 
tion is that there is a flux attached to each particle. In this way 
the phase associated to the original particle statistics is now 
generated by a mechanism based on the Aharonov-Bohm ef- 
fect. It can be seen that the uj function introduced by Fradkin 
1I29I1 (Eq. (I20» generates a vector potential A which is solu- 
tion of Eqjj26t. On the other hand, the lu function proposed 
by Wang |32| cannot be associated to any flux since the to- 
tal change of oj on any closed loop vanishes (or equivalently 
^ A.ds = 0). This result suggests that the vector potential can 
be eliminated by a gauge transformation. However, it is im- 
possible to find a gauge transformation which does not change 
the commutation relations of the particles. The advantage of 
the solution given by Wang fs^l is the straightforward gen- 
erahzation to higher spatial dimensions. Fradkin's approach 
(2d) i29il provides an alternative formulation of the original 
problem and has stimulated original methods to find relevant 
solutions in quantum Hall systems. The generalization of his 
approach to 3d isill is more involved since it requires an ex- 
tended Hilbert space and non-Abelian gauge transformations. 

Summarizing, we have shown that one can write down the 
generators of the canonical fermionic algebra as an operator 
function of hard-core bosons and then as a function of the 
generators of any Lie group having an irreducible represen- 
tation of dimension given by the dimension D of the local 
Hilbert space. If we are dealing with canonical fermions, 
the dimension D is equal to 2^f. However, we can obtain 
other values of D by imposing different constraints on the 
fermionic occupation numbers like in the bosonic case ll22il 
(see next Section). This is a common situation for strongly 
correlated problems where part of the Hilbert space is elim- 
inated when deriving effective low-energy theories. For in- 
stance, the Hilbert space of the effective low-energy theory 
for the Hubbard model [t-J model) has the constraint of no 
double-occupancy and therefore D is equal to 3. 



D. Anyons 

Similarly, one can extend this idea of transmutation of 
statistics to particles satisfying general equal-time anyonic 
canonical commutation relations defined by an angle 6. To 
this end we need to generalize the transmutators to any statis- 
tical angle < 9 < tt 



/3<a 



4 = exp[z-^c.(lj) 



nil 



(27) 



We have seen that there are many types of bosonic particles. 
For each type we will get a different type of anyon after sta- 
tistical transmutation. For instance, if we start from canonical 
bosons (uia = b\^b^J 

f[&i.,M = [^L4^]-0' (28) 



by transmutation of the statistics 

one gets anyons obeying commutation relations 



a. 



= 1 



(29) 



(30) 



and de/ormed commutation relations if both anyonic operators 
correspond to different sites or flavors. To write down those 
commutation relations we need first to define a particular or- 
dering for the combined site and flavor indices I33I1 . In this 

way, for(j,/3) > (i, a) 



. 



[ajp, flic 



(31) 



The requirement of an ordering for the indices comes from the 
fact that [A,B]g = — exp[i6][B, A]^0. So unless cxp[i6'] is 
a real number (i.e., bosons (9 = 0) or fermions (6 = tt)), we 
need to define a particular index ordering. It is in this order- 
ing that the intrinsic many-body character of the these parti- 
cles is encoded. In this simple example we can see that the 
Pauli exclusion principle and the phase associated to the ex- 
change of two particles with different indices (flavor or site), 
are two distinct and independent concepts. The exclusion be- 
havior is determined by the local commutation relations be- 
tween operators associated to the same site and flavor. Since 
these commutation relations are not changed by the transmu- 
tation, Eq. ( I29> , the exclusion properties are preserved when 
we change the statistics of the particles. In our example we 
can create any number of anyons in the same orbital and fla- 
vor because the local commutation relations are the same as 
the ones for canonical bosons. However, the deformed com- 
mutation relations of Eq. (13 U indicate that the result of ex- 
changing two anyons with different indices is the multiplica- 
tion by a phase factor cxp[i6']. It is this second aspect, not 
related to the exclusion properties, that decides whether the 
particles are bosons, fermions, or anyons. 

Now we will consider another example of transmutation 
from bosons to anyons. In this case, we will take the multifla- 
vored hard-core bosons defined by the commutation relations 
Eq. (II 1> . The particles obtained after the transmutation of the 
hard-core bosons will be called type I hard-core anyons 



(32) 



with 



7:^ = exp[z0 ^ , 



K? = exp[i-^cj(l,j) 711] 



(33) 



Like in the previous example, the local commutation relations 
are preserved {hja = aja^ja' = YHU "jc) 



0, 



(34) 
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In this particular case, since there is a hard-core condition 
aj^aj^ = 0, the operators also satisfy the following local an- 
ticommutation relations 



Again, we need to define an index ordering for the deformed 
commutation relations involving operators with different in- 
dices. For > {i,a) 



1 



(35) 




(40) 



Thus, the local anticommutation relations are also preserved 
under statistical transmutation. Clearly, Eqs. ( I35> are the lo- 
cal anticommutation relations for canonical fermions. This 
is not surprising since the multiflavored hard-core bosons de- 
fined by Eq. il 1> can be transmuted into canonical fermions 
(see Eq. il5\ ). For the commutation relations involving opera- 
tors with different indices, we have to define an index ordering 
like in the previous example. For (j, /3) > (i, a) 



0, 



(36) 



6 = corresponds to bosons which are hard-core in each 
flavor (see Eq. il and 9 = tt to canonical fermions, i.e.. 



.{9 = 0) 



tistical angle 9 one can put up to a single particle per mode 
and flavor (i.e, the particles are hard-core in each flavor and, 
therefore, satisfy the Pauli exclusion principle). It is evident in 
general that the transmutation of the statistics does not change 
the local commutation relations and therefore the exclusion 
rule. In a later section we will show how to generalize the 
exclusion principle to accommodate up to p particles per sin- 
gle quantum state. In particular, we will see that the para- 
particles, which obey this generalized exclusion principle, can 
be obtained from particles which are partial transmutations of 
canonical bosons (to obtain the parabosons) or fermions (to 
obtain the parafermions). 

The last example of anyonic particles corresponds to the 
transmutation of the hard-core bosons defined by the commu- 
tation relations Eq. @. This concept will be useful for the 
generalized JW transformations l22il that we introduce in the 
next section (these particles will be called type II hard-core 
anyons or simply JW particles) 



-J 



For any sta- 



(37) 



with 



I3<a 

K? = exp[i-^w(l,j) fii] 



(38) 



Notice that, formally, in all cases we could have considered 
different angles 9 instead of a single one in the expressions for 
the local, (e.g., 9i), and non-local, K? (e.g., 92), transmu- 
tators, although we do not see how physically relevant this 
general situation could be. 

It is important to stress at this point when the exchange 
statistics property of the particles (i.e., the property attached 
to the non-local part of the transmutator) becomes relevant in 
the description of a particular physical system. It is clear that 
whenever the system Hamiltonian does not permute particles 
on different sites i and j, then K? becomes a symmetry and 
the Hamiltonian is invariant under particle exchange statistics. 
This is the case in many Id problems (e.g., the XX Z model 
of Eq. or the Hubbard model) and in special 2d problems 
such as the U{1) gauge magnet that we expand on in section 

EXH 

Figure|2]summarizes the fundamental notions of languages 
and dictionaries connecting them, concepts developed in sec- 
tions |niand|in| 



Relation between modes: 

I Transmutation of Statistics 



Non-locality 



^ Locality j ^ 
f 

For each mode (local Hilbert space H of dimD ): 
A set of operators S acting irreducibly on H 



Language 



Set of operators S 
+ 

Statistics between modes 



Dictionary: 



( Language A ] (_ 



D 



B 



dim of the local Hilbert spaces are equal 



FIG. 2: What is a language, and when can two languages be con- 
nected? Summary of the main content of our fundamental theorem 
in conjunction with the concept of transmutation of statistics. 



IV. HILBERT SPACE DECOMPOSITION 



Since the local commutation relations are preserved, we have 
[aja:«jJ = [fljl^ajl] =0, ^^^^ 



We have seen in the previous section that the only condition 
to establish an isomorphism between two different languages 
is that the dimension of their local Hilbert spaces be equal (see 
Fig. 12). Notice, however, that the choice of local Hilbert space 
Hi depends upon the particular decomposition of the global 
Hilbert space H. When we work on a lattice, the most natural 
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decomposition is the one where Hi is the Hilbert space asso- 
ciated to each site (or mode). This is the decomposition that 
we have adopted in our previous examples. Another possible 
decomposition could be the one for which Hi corresponds to 
the Hilbert space of a bond (i.e., two sites instead of one). In 
general, there is no restriction in the way one can partition H. 

To be more concrete let us consider a 5*= 1/2 ladder system 
(see Fig.O. The global Hilbert space of the system can be de- 
compose into the direct product of the local Hilbert spaces on 
each rung. By doing so, the dimension of Hi is enlarged from 
D ~ 2, for the sites decomposition, to D ~ 4. According to 
our previous results, the change of D opens up the possibility 
of new mappings between languages. For instance, 13 = 4 is 
the dimension of the on site Hilbert space for canonical 5=1/2 
fermions on a lattice. Therefore, the 5=1/2 spin ladder can 
be mapped onto a system of electrons on a linear chain (see 
Fig- Ell. To be more explicit, a possible isomorphic mapping 
is (see Eq. ( I90t for a different mapping) 



= V2(5„ 



(41) 



Eq. 1411 can be inverted to obtain the expression for the spin 
operators as a function of the fermionic creation and annihila- 
tion operators 

'^j^ = Tstd - 2hil)cy<i + (1 - 2njT)A7cjJ , 



■Mi 



Q 



Q 



j+2 



Q 



j+3 



j+4 



■o— 














— •— 




i+1 


i+2 


i+3 



Bond j 



Site / 



FIG. 3: Mapping between an 5=1/2 ladder and a fermionic linear 
chain. In this case, bonds (rungs) are mapped onto sites. 



close the section showing that the notion of parastatistics in- 
troduced by Green l34ll is associated to partial transmutations 
of the canonical bosonic and fermionic languages. In section 
llXI we connect the gauge-field to spin languages through the 
quantum link model relation. In this way one envisions formal 
relations between all the different languages used to describe 
the properties of quantum systems. In the course of this sec- 
tion we will see that many well-known transformations are 
particular examples of applications of both the fundamental 
theorem and the transmutation of statistics. 



f 9+ - 

S!2 = 



T ^ 



(42) 



In this way we have mapped spin operators acting on a given 
bond onto fermionic operators acting on the corresponding 
site (see Fig. |3}- This dictionary will be exploited in section 
IVllto map ladder 5=1/2 magnets onto Id Hubbard-hke mod- 
els. 



V. BRIDGING THE LANGUAGES OF QUANTUM 
MECHANICS 

The purpose of this section is to illustrate with examples the 
algebraic framework developed in the previous three sections. 
In the first part we describe different SU{N) and SO{N) 
spin-particle mappings which are a direct consequence of 
our fundamental theorem. In these cases we are connecting 
bosonic languages and, therefore, it is not necessary to trans- 
mute the statistics. In contrast, in the second part we describe 
the generalization of the JW transformation where the trans- 
mutation of the statistics plays a fundamental role. In this 
case, we connect the different representations of 5t/(2)-spins 
to constrained fermions l^j. We devote the third part of this 
section to show how the fractional exclusion statistics alge- 
bras emerge from the present formalism in a natural way. We 



A. 517 (A^ ) spin-particle mappings 

1. Schwinger-Wigner bosons 

The connection between the Schwinger-Wigner bosons and 
the generators of su{N) in different irreducible representa- 
tions results from a simple application of the fundamental the- 
orem. We will consider a system of Schwinger-Wigner bosons 
with Nf ~ N different flavors (a G [l,iV]). These bosons 
obey canonical commutation relations 



and are characterized by the constraint 

N 

Q = l 



(43) 



(44) 



which sets the dimension of the local Hilbert space Hj, D = 
(^+//"^) . For = 2 we get 13 = M + 1, which means that 
we can get any value of D by varying M. We also know that 
the Lie algebra su{2) has one irreducible representation for 
each value of D. In this case D = 2S + 1, where 5 is the spin 
of the representation. According to the fundamental theorem. 
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if 25* + 1 = M + 1 we can write down the generators of su{2) 
{S^, , S^) in the S = M/2 representation as a polynomial 
function of the Schwinger-Wigner bosons and vice versa. The 
explicit form of these relations is the well-known connection 
between SU (2)-spins and the two-flavor Schwinger-Wigner 
bosons (a = 1, 2) 



5*? = 



St = 
SV = 



1 



(45) 



According to the fundamental theorem the inverse transfor- 
mation also exists. The creation and annihilation operators for 
Schwinger-Wigner bosons , Sj^ can be written as a polyno- 
mial function of S*?, SV , S?. However, since the Schwinger- 
Wigner bosons are not directly associated to the description 
of any physical system the inverse transformation is, appar- 
ently, not very useful. For this reason, the only application of 
the Schwinger-Wigner bosons is to provide a different frame- 
work for solving spin Hamiltonians. For instance, it is pos- 
sible to approximately solve the isotropic Heisenberg model 
by a mean-field approximation which preserves the SU{2) 
invariance of the model. This symmetry is violated in the 
usual spin-wave approximation where the spins are described 
by Holstein-Primakoff bosons. Through this simple example 
we see yet another important advantage of having different 
languages to describe the same problem: Different languages 
provide different frameworks to obtain approximate or exact 
solutions. Physics provides innumerable examples of applica- 
tions of this statement. It is not our intention to enumerate all 
of them in this paper, but to show the operative way to build 
the dictionaries connecting all of those different languages. 

To find the connection between the generators of the su{N) 
algebra and the Schwinger-Wigner bosons we just need to 
consider the general case of arbitrary number of flavors. To 
simplify the analysis we can start by considering the mini- 
mum value of A/, i.e.. A/ = 1 (i.e., the fundamental represen- 
tation). In this case the dimension of the local Hilbert space 
is D = (^) = N. Since TV is also the dimension of the fun- 
damental representation of su{N), the fundamental theorem 
states that the generators of su{N) can be written as a func- 
tion of the Schwinger-Wigner bosons. Again, this connection 
is very well-known and can be explicitly written as 



(46) 



Using Eqs. ( I43t and ( I44> we can verify that the 5"'^'s are the 
generators of the su{N) algebra in the fundamental represen- 
tation, i.e., they satisfy the following commutation relations 



(47) 



and operate in an A^-dimensional vector space. As in the 
N = 2 case, by considering larger values of M we obtain 
higher-order representations of su{N). The relation between 
the su{N) generators and the Schwinger-Wigner bosons is the 
one given by Eq. ( I46> independently of the representation. 



2. Fundamental irrep of su{N) and Jordan-Wigner particles 

The fundamental (quark) representation of su{N) can be 
mapped onto an algebra of constrained fermions (cj^, 

i 

Nf = N -1 flavors 



aj^l^w = tt)) or hard-core bosons (bj^ = '^la(^ ^ ^^^^ 





= a- aia 












Nf 
a = l 







where 1 < a, /3 < Nf runs over the set of particle flavors, 

^^''(j) (with < < 



and al„ = a]„n^ii(l- 
Nf) are the components of the SU {N)-spm (i.e., there are 
iV^ — 1 linear independent components). It is easy to verify 
that these are generators of an su{n) Lie algebra satisfying the 
commutation relations 



[5^^ (j)] = S^,,.S^'' (j) - (j) 

For instance, for = 3 we have Issll (a = 1,2) 



(49) 



SO) 



nn 



i 

nj2 



1 

3 



(50) 



We can immediately see that the 2x2 block matrix 
(1 < a, /? < 2) contains the generators of su{2). In general, 
from the commutation relations Eq. i49l . we can verify that 
if S^^d) are the generators of su{N), then are the 

generators of the subalgebra su{N — 1). This will be useful 
in section lVll 



3. Generalization to other irreps of su{N) 



Here we show that the hard-core bosons b^a can be con- 
nected to other irreducible representations of su{N). The di- 
mension of an su{N) representation of spin S = M/N is 



(N+M-l^ 
\ M 



The number of flavors is then N 



f 



/N+M-l\ 
M 



In this case, it is more convenient to adopt the following 
notation for the flavors: Each flavor a will be denoted by 
an array of N integer numbers {mf,m2, ■ ■ ■ satis- 
fying X]^=i "^CT ~ There are arrays sat- 
isfying this condition. However, there is a particular one. 



{rrii 



^0 



}, associated to the vacuum state. 



Nf 



a, 13=1 



(51) 
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Again, it is easy to verify that these are generators of an s 
Lie algebra satisfying the commutation relations Eq. (l49l 
This is just an example to show that the generators 
Lie algebra in a particular and arbitrary representation c 
mension D can be used to describe a system where D i 
dimension of the local Hilbert space. For each particular 
there will be representations which will be more approp 
(easier to handle) than others. For instance, this set of sr 
representations can be performed in a more natural way i 
Schwinger-Wigner bosons. 



B. SO (A'^) spin-particle mappings 

Another possible language related to the JW particles i 
S'0(A^)-spin language. By this name we mean the lang 
whose generators are the identity and the generators o 
so{N) Lie algebra (whose number is N{N — l)/2). It is 
to verify that the following set of operators 



M'^i) ^ {K?)^-a.^p--alKl 
X°°(j) = 0, 
obey the S'0(iV)-spin commutation relations 

-<^M-X'''''(j) + V-^'^'^O)], (54) 

with < ^, < TV and 1 < a, /3 < iV/ = iV - 1 The 
antisymmetric relation 7V1^''=— tM"^^ is immediately derived 
from the definition of M.^^" , Eq. \53\ . In this way, we have 
constructed another possible language connected to hard-core 
particles. Like in the previous case, the subset of operators 
M."^!^ is a set of generators for an so{N — 1) subalgebra of 
so{N). 

In the following, for simplicity, we will only consider 
the connection between 50 (A^) -spins and hard-core bosons 
0) with non-local transmutator K? = l. For 



instance, for N ^ i, the generators are (a = 1, 2) 
/ 

J J 





hi 



hi h; 



bib: 



bl b:. 



■}2 "i2 "}2"il "}l"i2 



jl"j2 "j2"}l 




(55) 



C. Generalized Jordan- Wigner Particles 

In 1928 Jordan and Wigner related the spin quantum me- 
chanical degree of freedom to spinless particles with fermion 



Classical spins to lattice gas or binary aJloy mappings 



Ising Model: (magnetism) 









• 


f J 

















{Si} e configuration; 5"; = ±1 



La,ttice Gas: (e.g., crystal melting) 



H 


>- 


H 


h- 











CoUection of "atoms" (no T) 
One atom at each site 



jV| jV, , jV| — > empty sites 



Binary AUoy: (e.g., ^- brass) 



Two types of "atoms" (no T) 
One atom at each site (• or •) 



FIG. 4: Classical version of the lordan-Wigner particles defined in 
l23l . Here the mapping is performed between the simple lattice of 
Ising spins Si — ±1 and a lattice gas of one type of atom or a binary 
alloy. Ncr, with a =t, |, is the number of Ising spins of type a while 
A*', is the number of hard-core atoms of flavor • with no kinetic 
energy T. 



statistics. This transformation involves the 5=1/2 irreducible 
representation of the Lie algebra su{2). We have recently ll23l 
generalized this one-to-one mapping to any irreducible spin 
representation of su{2), spatial dimension of the lattice and 
particle exchange statistics. From the physical viewpoint the 
JW particles are essentially hard-core independently of their 
flavor index (i.e., satisfy the Pauli exclusion principle with 
p ~ 2 (see below)). 

The generalized JW mappings constitute a kind of quantum 
version of the well-known classical spin-lattice-gas transfor- 
mations ll36ll (see Fig. 0}. In a classical lattice gas (or binary 
alloy) each site can be occupied by at most one atom, a hard- 
core condition equivalent to its quantum counterpart. On the 
other hand, since the spin system is classical there is no in- 
trinsic dynamics and the kinetic energy of the gas (T) must 
be included a posteriori in an ad-hoc fashion. This contrasts 
the quantum case where the dynamics of the JW particles is 
intrinsic. 

We have seen in section|^that the canonical fermions can 
be transformed into bosons which are hard-core in each fla- 
vor. We will consider now another type of fermions which 
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naturally emerge from the strong coupling limit of models 
for interacting electrons. If the short range component of 
the Coulomb repulsion is much larger than the kinetic energy, 
the repulsion can be effectively replaced by a constraint of no 
double occupancy. This perturbative approach is usually im- 
plemented by a canonical transformation, which leads to an 
effective Hamiltonian acting on the subspace of states with no 
double occupancy. The fermionic subalgebra used to describe 
this effective model is generated by the so-called constrained 
fermions. Therefore, the constrained fermions are obtained 
by imposing to the canonical fermions a local constraint of 
no more than one particle per orbital (or site). This constraint 
may be incorporated into the fermionic algebra by defining 
the following creation and annihilation operators for the con- 
strained fields 

4. = 4<.n(l - "j-) ' =11(1 - njO Cj, , (56) 



where J^^ is the set of flavors, with r/ = ^ , 1 depending upon 
the spin character of the irreducible representation. It is easy 
to check that the particles generated by this fermionic algebra 
satisfy the constraint of single occupancy, i.e., the eigenvalues 

of ~ 12a=i ^jcr either or 1. The most well-known 
context where these fermions appear in condensed matter 
physics is the strong coupling limit of the Hubbard model, 
which leads to the i- J Hamiltonian. While the Hubbard model 
is described in terms of spin- 1/2 canonical fermions (elec- 
trons), the fermionic language for the t-J model is generated 
by creation and annihilation operators of spin- 1/2 constrained 
fermions (ct, ct' =t, j) 



(57) 



In this section we will consider the general case of constrained 
fermions with Nf = 25* different flavors. In the context of our 
previous example, this generalization can be interpreted as the 
natural language for a t-J model with more than one orbital 
per site j or larger-spin fermions. The set of commutation 
relations for only two different flavors, Eq. i57\ . is generalized 
in the following way 



{q,,Cj,,} = {cl,cl,} = 0, 



1 + Jijo- - nj 



if CT = ct', 
if CT 7^ ct'. 



(58) 



Uj with number 



Notice that Yir'eJ^ ^ '^J^) ^ 1 + '^jo- 
operators satisfying fi-^^n-^^i = ^o-o-'^-jo-- 

We may now ask what is the bosonic language obtained 
from the constrained fermions when the statistics is trans- 
muted, i.e.. 



v. 



}a J 



(59) 



(We do not need to include the local transmutator since the 
constraint does not allow to have more than one particle per 



site (or orbital).) The answer to the question becomes clear if 
we compare the set of commutation relations for the hard-core 
bosons, Eq. (|6}, to the set of anticommutation relations for the 
constrained fermions, Eq. (I58> . All antic ommutators which 
are not of the form {cj^jcj^,} are mapped onto the corre- 
sponding commutators for hard-core bosons, Eq. (|6j. This re- 
sults from the transmutator KJ . What is the effect of the trans- 
mutation on commutators containing bilinear forms in the cre- 
ation and annihilation operators? First, it is easy to check 

that in this case the commutator [^L^i/S'^] ~ '^ij'^07^L °f 
bosonic operators is mapped onto the corresponding commu- 
tator of fermionic operators 



(60) 



Finally, to complete this correspondence between commuta- 
tion and anticommutation relations we have to consider the 
effect of the transmutation on the commutators b\^] and 

the anticommutators {Cjg.,Cjg.}. It is easy to check that the 
transmutation does not change these particular products, i.e.. 



(61) 



We can see that the anticommutator {6;^, b\^} is a linear com- 
bination of generators of the Lie algebra defined by the set 
of commutation relations, Eq. This is so since the set 
of generators {/, b\^, b^^^b\j}-^p} belongs to the fundamental 
representation of the Lie algebra. 

Therefore, in this case, the generators of the Lie algebra 
associated to the bosonic language are transmuted into gen- 
erators of a Lie superalgebra associated to the constrained 
fermions. The Lie product in the Lie algebra, Eq. is 
turned into the graded Lie product in the superalgebra. The 
generators of the bosonic Lie algebra can be separated into 
two subsets which are mapped onto the odd and the even gen- 
erators of the fermionic Lie superalgebra. In this particular 
case, the Lie algebra associated to the commutation relations, 
Eq. (|6j, is ii(l) su{N). The creation and annihilation op- 
erators for the hard core bosons, b\^ and 6;^ are mapped onto 
the odd generators of the Lie superalgebra, and q^, while 
the identity and the biUnear forms, / and b\pb-^^ are mapped 
onto the even generators / and Cj^Cj^. The bosons which are 

hard-core in each flavor, 6;^, and the canonical fermions, C;^, 
provide another example of transmutation of a Lie algebra into 
a Lie superalgebra. In this case, the generators of the Lie al- 
gebra u(l) su{2) are transmuted into the generators of 
a Heisenberg Lie superalgebra. 

In this way, through the transmutation of statistics we have 
established a direct connection between the multiflavored con- 
strained fermions and the multiflavored hard-core bosons. Us- 
ing this connection we can see that the MM transformation 0| 
can be obtained from the JW transformation [3J by a simple 
transmutation of the statistics. 

In the course of demonstrating the Corollary I, we have 
shown that for each class of languages characterized by the di- 
mension D of the local Hilbert space, there is one whose gen- 
erators are also generators of the Lie algebra u{l) su{2) 
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(D = 2S + 1). In addition, we have seen that there is another 
language in the class whose generators are the identity and 
the creation and annihilation operators for multiflavored hard- 
core bosons {Nf = D — 1). Then, according to the fundamen- 
tal theorem, the three components {Sf, , S?} of a spin S 
can be written as a polynomial function of {/i, bl^,b^^} with 
1 < a < Nf and vice versa. This is the generalization of the 
MM Qftransformation to any irreducible spin- 5 representa- 
tion of su{2). Adding now this result to the mapping already 
established between the multiflavored hard-core bosons and 
the constrained fermions, we can conclude that it is also pos- 
sible to write {S?, , S^} as a polynomial function of the 
identity and the creation and annihilation operators of multi- 
flavored constrained fermions and vice versa. The later is the 
generalization of the JW transformation to any spin S j23l . 
On the other hand, the generalized JW spin-fermion mapping 
can be easily extended to include a spin-anyon mapping sim- 
ply by using the anyonic particles generated by at and a- 



(see Eq. ( BTt l 
given bv Il22ll3" 



MM- The explicit form of this generalization is 
Half-odd integer spin S (a e J^i = {-S + 1, . . . , 5^}): 



4. 



2 

-S+ ^ (5 + ct) nja , 
{K?)^L3 (S+) V? , 



where 



T-Sf 



cr-1 



rey^i T= — S 



Integer spin S {a e Ti = {-5, . . . , -1, 1, . . . , S}): 



(Te:Fi 

\[S,) VI ifa<0, 
where Vj =\[ ^1 = \[ 77-1 



and 77cr = \/{S - a){S + fi + 1) (see Fig. |5}. [A bar in a 
subindex means the negative of that number (e.g., a = —a).} 



Flavors 

S — 

S-\ 

S'-2 



5^ 

s 
s-\ 

S-2 



Flavors 

1 

— S 



-1 

-S + 2 
-S + \ 
-S 



O 



-1 



-s+2 
-S + 1 
-S 



S-1 
s-2 



-S+2 
-S + 1 



Integer S 



o 

Half- odd Integer S 



FIG. 5: Constrained any on states per site for integer and half-odd in- 
teger spin S. In both cases there are 25* flavors and the corresponding 
25 + 1 values of 5*^ are shown in the middle column. One degree 
of freedom is assigned to the anyon vacuum (circle) whose relative 
position depends upon the spin being integer or half-odd integer. 



The total number of flavors is Nf = 2S, and the 5=1/2 case 
simply reduces to the traditional JW transformation. These 
mappings enforce the condition on the Casimir operator = 



8(8 + 1). 



D. Fractional Exclusion Statistics Algebras 

What microscopic properties of the fundamental particle 
constituents determine the thermodynamic behavior of mat- 
ter? It has been known since the early days of quantum me- 
chanics that particles in nature are either bosons or fermions in 
the three-dimensional space. It was immediately realized that 
for a given interaction among particles the properties of matter 
were completely different depending upon the constituents be- 
ing bosons or fermions (for example, normal Coulomb matter 
is stable thanks to the fermionic nature of its particles). More- 
over, it has been known for a long time that there are other 
theoretical possibilities that interpolate between fermions and 
bosons, i.e., cases where p > 1 identical particles occupy one 
and the same state. Quantum statistics was the name given to 
this phenomenon clearly nonexistent in the classical physics 
description of matter 

In the seventies it was realized that two space dimensions 
allowed for other exotic possibilities with particles having 
fractional spin and statistics f39fl. Later on, this concept 
proved to be not merely an intellectual exercise since parti- 
cles obeying such quantum statistics were experimentally re- 
alized in the context of the quantum Hall effect (QHE). It was 
Haldane j40ll . to our knowledge, the first to realize that the 
Pauli exclusion principle and exchange statistics ideas were 
in principle independent concepts. In this way, motivated by 
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the properties of quasiparticles in the fractional QHE, he gen- 
eralized the Pauli exclusion principle introducing the concept 
of fractional exclusion statistics. 

Basically, two physical features characterize the quantum 
statistics of particles: One is related to their indistinguishabil- 
ity and is characterized by the property that when two iden- 
tical particles are exchanged the total wave function acquires 
a phase factor (exchange or permutation statistics). Another 
refers to the ability to accommodate p particles in the same 
single-particle quantum state (exclusion statistics) EH. The 
first concept depends upon the space dimensionality of the 
system while the second one is independent and, therefore, 
unrelated to the notion of anyonic fractional statistics which is 
applicable exclusively in two space dimensions. For instance, 
from the generalized Pauli principle viewpoint there is no dis- 
tinction between particles obeying the algebra of Eqs. (I39> and 
(I40> for any statistical angle 6, i.e., all of them correspond to 
p — I. However, the exchange of two particles a does depend 
upon 9. 

In this section we consider the problem of formulating the 
algebra satisfied by these generalized Pauli-exclusion parti- 
cles and its connection to our fundamental theorem. Our def- 
inition is a possible second-quantized version of Haldane's 
definition introduced at the first quantization level. 

To simplify matters we consider as an example a single fla- 
vor (single site) fractional exclusion statistics algebra 



9U 



1 -i^ , 



= 0, 



7t 



(62) 



where the operator F (F^ = i^) is a polynomial function 
of and g , such that g^F = 0. The F term represents a 
sort of deformation of the canonical boson algebra (deformed 
Heisenberg algebras play an important role in the theory of 
representations of quantum groups). For example, F can be 
= ia^Yia y = 11^=0(5^5 - j) and p is an arbitrary 
integer with the condition p > 0. For g-particles with this F 
operator it is very easy to prove that they satisfy the nilpotency 
condition 



(.9 



= 



(63) 



which means that one can put up to p particles on each mode. 
A {p+ 1) -dimensional matrix representation of this algebra is 



/DIG 
V2 
\/3 




\0 






Vp 
0/ 



(64) 



Certainly, one can build infinitely many different algebras 
leading to different types of fractional exclusion statistics par- 
ticles with the same value of p. For example, consider the 
su{2) algebra generated by the operators {gl,gi,n?} (with 



K) 

i.e.,p = 2) 



I?) in the S=l (triplet) representation ({gj 



t^3 



= K,.9]] =0, 
[.9i, ] = - nf) , [n?, gj] = gJ 



0, 



(65) 



It is straightforward to realize how these particles relate to the 
usual spin generators of su{2) 





= V2gl, 


Si 









(66) 

which can be seen as a possible generalization of the MM 01 
transformation for S=l. 

For example, if we consider the Id XY model in the S=l 
spin representation 

i/xy = jY^iSfSf^, + SfSf^,} , (67) 
j 

it is simply related to the g-particle Hamiltonian 

H,y = ^I](5l.9j+i+,9j+i5j) 
j 

= 2J^coskgJ,5k, (68) 

k 

where §1^ = (l/\/]V7) exp(ik • j) gj is the Fourier- 
transformed operator 

E. Green's Parastatistics 

Since the very beginnings of quantum mechanics peo- 
ple tried to understand the depth and consequences of the 
symmetrization postulate which asserts that physical states 
of identical particles must be either symmetric (Bose statis- 
tics) or antisymmetric (Fermi statistics) under permutations. 
In 1953 H.S. Green |34] considered the possibility of hav- 
ing non-identical but dynamically similar particles satisfying 
what is now known as parastatistics. The introduction of 
parastatistics allowed some hope to represent all particle fields 
in terms of a fundamental spinor field and, of course, canon- 
ical fermions and bosons represented particular cases (paras- 
tatistics of order p = 1). 

Soon, Greenberg j43l suggested that quarks satisfy paras- 
tatistics of order 3 and Green f43l reformulated the neutrino 
theory of light identifying the neutrinos with parafermions of 
order 2. It is not our intention to discuss the physical aspects 
of these assertions but to present parastatistics as another op- 
erator language which allows representation of Lie algebras 
and to relate these particles to the ones already introduced in 
previous Sections. 

Following Green l34l l43ll . let us introduce the auxil- 
iary modes (upper (lower) sign will be used to define the 
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parafermions (parabosons)) 



(69) 



labeled by the Green indices a = 1, - ■ ■ ,p and j = 
1, • • • ,N and where df = {dfY, together with the condition 
I vacuum) = , V(j,a). For a 7^ /3 the auxiliary modes 
satisfy non-standard relations, i.e.. 



K,d?h = [d?J^]^ = 



(70) 



Parafermion and paraboson creation and annihilation opera- 
tors are defined in terms of the auxiliary modes by 



and consequently satisfy the commutation relations 



[K,dj]^,d,]_=0. 



(71) 



(72) 



It easy to verify that one can define a set of commutative num- 
ber operators (71? = d^d^) by 



nf = -{[dld,]^±p) 
which satisfy the commutation rules 



p 



a 



(73) 



(74) 



Note that in the parafermionic case (n?) = nj", which im- 
plies that nf has eigenvalues ranging from to p. Moreover, 



(4 



p\ 



114 



(75) 



such that {dl)P+'^ = 0, which in a sense also generalizes the 
Pauli exclusion principle. It is clear from these definitions that 
canonical fermions (bosons) are parafermions (parabosons) of 
order p = 1. 

To connect these para-particles to the ones already de- 
scribed in previous sections, one has to realize that (total or 
partial) local and non-local transmutations are necessary. It is 
easy to see that the parafermionic auxiliary modes are 



d? 



j<i I3^a 



(76) 



where a partial non-local transmutation connects to the 

canonical multiflavored fermions cj^, Eq. ilQi . Similarly, the 
parabosonic auxiliary fields are obtained as 

df = bl exp[*7r(^ E "J/' + E "i/^)] ' (77) 

j<i I3^a I3<a 



in terms of canonical multiflavored bosons fej^, Eq. ( I28> . 

One can represent all Lie algebras in terms of parafermion 
or paraboson creation and annihilation operators. Several au- 
thors contributed to this observation in the early sixties and 
seventies 13 ■ In the following we just present as an example 
the isomorphism between the algebra of N parafermions and 
the proper orthogonal complex Lie algebra Bm (i.e., the com- 
plexification of so(2iV + 1 )) of dimension TV (27V + 1 ). The set 
of operators antisymmetric in their indices (m, n e [1, 2A^]) 

^m^n) , (m < n) 
LmO = ^^Om = i (78) 

^00 = , 



where 



d\ + d. d\- d. 



2i 



(79) 



represents a basis which span the Lie algebra isomorphic to 

\L ^ Li,u''\ — i{^—8^'yL^i/' -\- S Iff L 

+ (^aiiz-^mV - (^AiL/'i/iv) • (80) 

VI. EQUIVALENT CLASSES OF MODELS 

In the previous section we have shown how the different 
languages of quantum mechanics are connected. How can 
we exploit these dictionaries to get a better understanding of 
different physical phenomena? There are many different an- 
swers to this question which we will develop in the next sec- 
tions. In particular, the present section is devoted to show 
the equivalence between models which, in principle, describe 
completely different systems 14^ . In other words, we show 
that models describing different physical problems are asso- 
ciated to the same Hamiltonian written in different languages. 
In this way, these dictionaries appear as an essential tool to 
connect distinct fields of physics. Preliminary applications of 
these concepts are given in Refs. |2^ l46ll . In Ref. l23l . we 
showed explicit connections between spin S=l and t-J-like 
models. In Ref. |46] we demonstrated the equivalence be- 
tween SU{3) ferromagnetism, coexistence of magnetism and 
Bose-Einstein condensation in an interacting bosonic gas, and 
coexistence of a spin-nematic phase and ferromagnetism for a 
spin one system. The choice of examples in the present sec- 
tion [Vl] is by no means exhaustive, it is simply illustrative. 
Table |I] summarizes some of the most celebrated mappings 
known today. 

Another interesting aspect of these mappings is the pos- 
sibility of unveiling hidden symmetries which lead to exact 
or quasi-exact solutions. By exact solvability we mean that 
the full spectral problem is reduced to an algebraic proce- 
dure, while problems which are quasi-exactly solvable are 
those which admit only a partial determination of the spec- 
trum by an algebraic procedure (see Fig. |6}. Sometimes there 
are particular languages which allow us to recognize invari- 
ant subspaces of our Hamiltonian. The restricted action the 
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[/(I) gauge magnet 


F Strings 


c- 


E 


2 


Shastiy-Sutherland (5=1/2) 


t-t'-V-V 


c- 


GS 


Any 


BB 5 = 1 ((^ = 57r/4) 


t-J-V-fi 


c- 


QE 


Any 


Anisotropic 5 = 1 Heisenberg 


B (5=1/2) Hubbard 


P- 


U 



TABLE I: Equivalence between spin models (model A) and parti- 
cle models (model B), equivalence that can be simply established 
by using our fundamental theorem. The meaning of the different 
terms is clarified in the text (BB: bilinear-biquadratic, F: fermionic, 
B: bosonic). Last column describes the type of equivalence (C: com- 
plete Hilbert space mapping, P: partial Hilbert space mapping) and 
type of solutions (E: exact, E (BA): exact Bethe-ansatz, GS: only 
ground state, QE: quasi-exact, U: unsolvable). 



Hamiltonian in an invariant subspace can always be described 
by a class of languages which are more basic than the original 
one. By more basic we mean a lower dimension of the local 
Hilbert space, or in other words, languages describing sys- 
tems with less degrees of freedom. There are particular cases 
where the simplification is even deeper, since by using a more 
basic language we also recognize hidden symmetries for the 
restricted action of the Hamiltonian in the invariant subspace. 
The quasi-exact solution of the Id t-J^ |0] model, is an exam- 
ple of successful application of this concept. This quasi-exact 
solution led to the exact quantum phase diagram and charge 
excitations of the t-J, model M. 



A. 5(7(2) Heisenberg Magnets 

The XYZ 5'= 1/2 spin model {D = 2) is the most popular 
example of a family of Hamiltonians which is Bethe-ansatz 
solvable in d= 1 L47J 



(iJ> 



T QV QV 
Jy O; Oj 



Jz S-j Sj . 



(81) 



Using the MM 0] transformation of Eq. (|8}, the XYZ 
model is mapped onto a gas of interacting hard-core bosons 
with density-density interactions and particle non-conserving 



Complete Hilbert space 
Model A Model B 




Quantum isotropic XY 
Bilinear-biquadratic 5=1 



Partial Hilbert space 
Model A Model B 







— t-J, chain 

— Bosonic (5=1/2) Hubbard 



FIG. 6: Schematics of the two types of mappings that one can per- 
form: Complete and partial Hilbert space mappings. The figure also 
shows examples illustrating each class. 



terms 



(iJ> 



(iJ) 



(82) 



(ij) 



where t ^ {Jj, + Jy)/A, A = (J^, - Jy)/4, and V = J^. 
Note that in d=\ the non-local transmutator is a symmetry 
and, therefore, we could have used spinless fermions instead 
of hard-core bosons in Eq. ( I82l i. 

Another family of spin Hamiltonians {D = 3) which has 
been extensively studied is the bilinear-biquadratic (BB) S=\ 
Heisenberg model ( J > 0) 



= J\/2^[cos0Si • Sj +sin0(Si • Sj 

(iJ) 



(83) 



Using the S=\ generalization of the JW transformation for 
hard-core bosons ||22 



V2(6j,+6],) 

nji - nji (84) 



Si 



(see section rVCt . i/^ is mapped onto a i-J-like model for 
5=1/2 hard-core bosons including particle non-conserving 
terms 

^0 = E '^Ih^ + ^ "^Li + H-c.) + 4A ^ 5? . s? 

(i.D.'T (ij) 
(ij) (ij) j 



where s; 



cj — ^bj^cTafsb^p (cr denoting Pauli matrices), t 

smd), V 
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J, /i = zJ\/2sin0, and z is the coordination of the lattice. 
Again, for d=l the exchange statistics of the particles is irrel- 
evant and one could have used S=l/2 constrained fermions 
instead of hard-core bosons. From the known solutions in 
the spin model one can immediately recognize the solvable 
cases in the particle model Eg . J85> . These are: (f> — -k/A 
{d=\) Lai-Sutherland (LS) H, = 77r/4 (d=l) Takhtajan- 
Babujian (TB) tan(/) = 1/3 (d=l) Affleck-Kennedy- 
Lieb-Tasaki (AKLT) Q, 37r/2 (d=nKlumper (K) isill . 
and = 57r/4 (any d) our work in Ref. 

Of particular current interest are the 5=1/2 Heisenberg 
models on ladders. The simplest case corresponds to having 
only nearest-neighbor magnetic interactions (see Fig.Qi) 




i+1 



i+2 



i+3 



i+4 



Using the mapping given in Eqs. (I42> which connects 5=1/2 
spins with canonical fermions, -ff^gj^ can be rewritten as a 
Hubbard-like model on a linear chain (up to an irrelevant con- 
stant). For instance, if Ai = 1 and A2 = we get 

-[/E"jT^ji, (87) 
j 

where U ^ Ji and t = J-z- This is a correlated Hubbard 
model with zero two-body hopping terms plus a supercon- 
ducting term. In the absence of the superconducting term, 
this model has been exactly solved by Arrachea and Aligia 
ll52il . To eliminate the superconducting terms the original spin 
model has to be modified in the following way 

j 

0.1^=1,2 

where Pj = E,.=i,2('5jl + '^j+i..) "JT ^ ^ji + "j+iT - 
fij+il = 2{s?+s?_^_^). H^^^f^ for Ai = 1 and A2 = is equiv- 
alent to the Hubbard model with correlated hopping solved in 
Ref. ll52il . If the sites 1 and 2 are interchanged in one sublat- 
tice (see Fig.0>), the sign of the three-body hoppings changes 
from positive to negative. 

Let us consider now the Heisenberg spin-ladder model il- 
lustrated in Fig.|8l 

Hss - Ji E(^i^i^j2 + 5?2 + SISD 
j 

+J-2^{^2Sj^Sj^i^, + Sj^sf_^i^, + sy^sy_^_-^^^,). (88) 



FIG. 7: A spin S=l/2 Heisenberg ladder mapped onto a Hubbard 
chain model with correlated hopping. Models a and b are related 
through a gauge transformation. 



The isotropic limit (Ai = A2 = 1) of this model has been 
considered by Sutherland who established that for Ji > 
2 J2 the ground state is a product state of bond singlets. Again, 
using Eqs. ( I42> . Hss can be rewritten as a Id Hubbard-like 
model (up to an irrelevant constant) 



Hss = ^E 4-4+1^(1 - ^J-)(l - + 

+ ^E 4'T^j + l<T(l - njff)(l - flj+la) + H.C. 



where U ~ JiAi, t = J2, Jz = 4J2A2, and fj, = 
(1 — Ai)Ji/2. This model is quasi-exactly solvable in the 
isotropic case Ai = A2 = 1. In addition, if we can- 
cel the superconducting terms, the resulting model is also 
quasi-exactly solvable because there are invariant subspaces 
for which the action of the Hamiltonian can be mapped onto 
a t-Jz chain model 0| (remember that a constrained 5=1/2 
fermion is related to the canonical ones through the relation 

For correlated hopping model Hamiltonians, the previous 
mapping, Eq. ( I42> . works well. However, for standard (non- 
correlated) Hamiltonians it does not. Let us introduce a new 
transformation that is expected to work in the latter case 



(90) 







<% 


= ^jT - 1 


^^2 




cz 
K'='j2 
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i+1 



i+2 



i+3 



i+4 



FIG. 8: Sutherland 5=1/2 model and its mapping onto a Hubbard- 
like chain model. 



where the nonlocal operator K^a is now defined as 

i^jl = exp[i7r(^ fill +^ nit)] , 
1 Kj 



(91) 



Kj 



This transformation was originally introduced by Mattis and 
Nam 1.541 to solve the following Hubbard-like model 



which after the spin-particle transformation, Eq. J90L be- 
comes (i^ = 1,2) 

Af-l N 

Hmn = X/ ^i'^^i+ii' + ^ X/ '^ji'^j2- (93) 

Using the same transformation, the one-dimensional Hubbard 
Hamiltonian 

N-l ^11 



j=i 



2" 2' 
(94) 



can be transformed into 

N-l N 
-^Hubb = 2t ^ ('5'jV'5j + i^ + Sj'^Sj'|_i^) + U ^ 'Sjl<Sj?2, 

j=l 

(95) 

which represents a two-leg ladder made out of two XF-chains 
coupled by an Ising interaction. 

Other interesting examples of quasi-exactly solvable mod- 
els in 2d of relevance for strongly correlated matter will be 
presented in a separate publication Issll . 

B. Single-band fermionic Hubbard Model 

The Hubbard model is the most popular model of a strongly 
interacting system in condensed matter physics. It contains a 



kinetic energy term represented by a hopping integral t plus a 
local on-site Coulomb repulsion U. The single-band Hubbard 
Hamiltonian is {a I) 



(96) 

The dimension of the local Hilbert space is 13 = 4 
since per lattice site j we can only have the states: 
{ |0), cj^ |0), cjjO), Cj^cjjO)}, and |0) is the vacuum state. Re- 
placing the c-operators by the expression for transmutation of 
statistics given in Eq. ( I15l l we get 



(97) 



where 7T| = exp[z7rriij, Tj^ = 1, /ij = cxp[i a(l, j) n\] 
with fij = fiji + njj^, and we have used the result 



(98) 



The (statistical gauge field) vector potential associated to this 
transmutation of statistics is 

1 

Using this expression for the vector potential we can rewrite 
^Hubb in the following way 

i?Hubb = t Y (^LcxpM'.(i)] V +H.c) 



(100) 



In this way we see that the Hubbard Hamiltonian can be writ- 
ten in a bosonic representation but with interactions which are 
non-local when d > 1. When = 1 the vector potential acts 
as a correlated hopping term and the interactions become local 

^^Hubb = ^Et ^iT(l - 2ni+u)6i+n + blil ~ 2niT)&i+u 

i 

+ H.c] + C/^(niT-i)(nu-i). (101) 

The 5J7(2)-spin 5=3/2 is another physical object for which 
the dimension of the local Hilbert space is = 4. From 
our fundamental theorem it is possible to write down these 
constrained bosons in terms of spins 5=3/2. In this way, 
we can find a 5=3/2 representation for the Hubbard model. A 
possible mapping between these bosons and a spin 5=3/2 is 



5. 



3 1 



Sr = ^/3 6j, + 2 5],6j,, 
3.. 



5f 
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This mapping can be inverted to get the bosonic operators as 
a function of the spin operators 



1 



2\/3 
1 



+ \2 



^/3 



(103) 



Similarly, the expression for the spin 5=3/2 operators in 
terms of canonical fermions is 



s; 



4, 



2("JT 



-jVjT ' 
o("jT 



1) 



(104) 



Again, we can write down the fermionic operators in terms of 
the spin operators 



1 



(105) 



Using these expressions we can write down the Hubbard 
model in terms of 5=3/2 spins as 

(i-j) 



2 

-E 

4 ^ 



H.c] 



(106) 



The 5=3/2 representation of the Id Hubbard model is 

'^[5+(5| + i)(5?,i 



Id 

Hubb 



3 
J 



(107) 



Since the Id Hubbard model is Bethe ansatz solvable Is^l . we 
are providing a new spin 5=3/2 Hamiltonian which is also 
Bethe ansatz solvable through the isomorphic mapping be- 
tween the two languages (fermions and spins). 

We can also write down the Id Hubbard model using the 
hierarchical 51/(4) language. To do so, we fist need to find the 
dictionary connecting the components of an 5[/(4)-spin in the 
fundamental representation with the creation and annihilation 
operators of the particle language 



(l-nji)fe]^ 5"(j) 



(108) 



where = (1 - - nj^) - i, - ^jT(l - 

^ji)-i'522(j) = njj^(l-nj|)-i,and533(j) = nj|7ij^-i. 
For a general SU{N) group there are actually two kinds of 
spinors: upper and lower. The upper spinors transform ac- 
cording to the conjugate representation. For the particular 
case of SU{2), the conjugate representation is equivalent to 
the original one, i.e., the conjugation is equivalent to a sim- 
ilarity transformation. However, in general for N > 2 the 
conjugate representation is not equivalent to the original one. 
Consequently, the SU{N > 2) ferromagnetic and antifer- 
romagnetic Heisenberg Hamiltonians are essentially different 
operators (in the case of SU{2), they just differ by an overall 
sign) iSTil . In the same way we wrote in Eq. jl08t the genera- 
tors of su(4) in the fundamental representation, we can write 
down the corresponding expressions for the generators in the 
conjugate representation 



5(j) = 



("JT - l)^ji 
("ji - l)^jT 



-522 (j) 

(109) 

i, and533(j) = 



where = njinji 

(j) = njt(l-nji) 

"ji) - 3- 

The expression of the Id Hubbard model in terms of the 
SU{4) hierarchical language is ( J^^, J^^, J^^ = J^^) 



TTld 

^Hubb- 



1) 



U 



(110) 



where the non-zero magnetic interactions are: Jqi = J02 = 

—Ji3 — —^23 = i and Jqx -J02 = ^-^13 = ^"^23 = 
—t. We can see from Eq. ( II 10> that in this representation the 
U term plays the role of a magnetic field. The case ,/^^ = 
corresponds to the Hubbard-like model which was exactly 
solved by Arrachea and Aligia Is^l . 



C. Single-band bosonic Hubbard Model 

We consider now a model of hard-core bosons including 
two degenerate orbitals a = {1,2} per site 

(ij),a,a' j 

+ y5](ni-l)(nj-l)-ME^J' (111) 

(iJ> j 

where hja = bjcfija ^^'^ "j = "ji + "j2- The first term 
represents a hopping connecting pair of orbitals belonging to 
nearest-neighbor sites The second is a local Coulomb 

repulsion term between particles occupying different orbitals 
of the same site. The third term corresponds to a nearest- 
neighbor density-density repulsion V. The dimension of the 
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local Hilbert space is D ~ A since per lattice site j we can 
only have the states: {\0) ,bl^\0) ,1^0) ^bj^byo)}, and |0) is 
the vacuum state. 

As mentioned in section IIIl Bl the hard-core boson opera- 
tors in each site generate an su{2) algebra: {^iq., ^jQ, '^ia — 
1/2}. The associated representation is the fundamental one, 
i.e., 5=1/2. Since there are two orbitals per site, the local al- 
gebra associated in each site is the direct sum of two su{2) al- 
gebras, and the representation is the direct sum of both 5=1/2 
representations. It is well known that by reducing the direct 
sum we get two irreducible representations: the singlet and 
the triplet representations. The singlet is associated to the an- 
tisymmetric (under the permutation of both orbitals) state 



1 



\^iA)^^ibl-b],m^bl\0) 



(112) 



The three remaining states, which belong to the triplet repre- 
sentation, are symmetric. These states can be generated by the 
creation operator 



1 



V2 



(113) 



By applying this operator to the vacuum state |0) one can get 
a particular basis for the triplet representation 



{|0>,6ls|0),6ls6y0)} 



(114) 



The singlet plus the triplet states generate another basis for the 
local Hilbert space of H with L» = 4 (nj = bj/J)jA + ^js^js)- 
If we apply the MM 0] transformation to these hard-core 
bosons, we get a spin 5=1/2 for each orbital. After that 
transformation it becomes clear that ji/'iA) is the singlet state, 
while the other three (Eq. il 14» are the triplet states with 
52 = —1,0,1. Note that the local algebra satisfied by 6!^ 
and bl^ is not the same as the one satisfied by bj^ and fe;,. In 
particular, P-J)}^ = bj^blj^ = 0. 

Going back to our model, Eq. il 1 1> . we immediately no- 
tice that the local singlet state I^/'ia) is invariant under the ap- 
plication of H. In other words, if there is one particle in an 
antisymmetric state at site i, that local state will be conserved. 
There is a local U{1) symmetry associated to this conserva- 
tion 



(115) 



which leaves H invariant. Consequently, one can identify in- 
variant subspaces of H. These subspaces are classified ac- 
cording to the set of sites A/a which are in a singlet state. We 
can now take advantage of the invariance of these subpaces, 
and project H onto each of them to reduce the effective num- 
ber of degrees of freedom. The expression for H (up to a 



constant) restricted to the subspace associated to the set A/a is 

(ijy J 



+ Z^'^jsl'^js - 1) , 



(116) 



where and j means sites i,j ^A/Aandnjs = 1 — [^jsi^jsl- 
Having a singlet at a given site j is equivalent to take the site 
out of the lattice since that local singlet is frozen at j (Pauli 
blocking). The action of H on each invariant subspace can be 
described by the operators fojg and b^^, which are generators 
of an su{2) algebra in the 5=1 (triplet) representation. This is 
easily seen using the mapping of Eq. ( I66t after replacing the 
^j's by the fejg's operators acting only on sites j ^ A/a 

5+ = V2bl, 
Sr = ^/25js, 

5? = njs-l. (117) 

A natural consequence of this transformation is the possi- 
bility of mapping iJvA onto an 5=1 spin Hamiltonian. This is 
done by inverting Eq. ([^2} 



njs = 5/ + 1, (118) 
and replacing these expressions into Eq. (II 161 to get 



1 



(ijy (ijy 



(119) 



This is an anisotropic Heisenberg model ( = V and = 
Jy = —2t) with a magnetic field = U /2 — ji applied 
along the z axis. In addition, there is an easy-plane single- 
ion anisotropy. Therefore, each invariant subspace of H gives 
rise to an anisotropic Heisenberg model Hf^^^ acting on a par- 
tially depleted lattice (i.e., the sites in A/a are removed from 
the lattice). From the variational principle it becomes evident 
that the lowest energy subspace is the one containing no local 
singlet states, i.e., Aa is the empty set (the original lattice is 
not depleted). 

We know from our fundamental theorem that there are 
many other possible languages that may be used to describe 
the present system. In particular, since (b\^Y' = it is clear 
that these modes satisfy a generalized Pauli exclusion princi- 
ple with p — 2 and, therefore, one can establish a mapping to 
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g-particles which satisfy the algebra written in Eq. i62\ . The 
explicit form for this mapping is 



1 



~V2 



njs 



■V2 



(120) 



whe re n? = gjgy Again, by replacing these expressions into 
Eq. ( II 16l l. we can write down in terms of the g-particles 



as 



H^^ = -Y.^9lgj + glgi)ihi + h2 + h3) 

(ijy 



(iJ> 



(121) 



where hi = ti{n? + — 2){n? + — 3), /i2 = ^t2{n? + 
n? - 1 ) (nf + n? - 3), and = (nf + n? - 1 ) (nf + n? - 2) 

with ti — t,t2 = V2t, and = t/2. The correlated hopping 
terms values are such that the matrix elements of the three 
possible hopping processes, illustrated in Fig.|9] are the same. 



between modes in Fourier space, for instance. A very well- 
know example is the one introduced by P.W. Anderson Is^l 
for the case of the Bardeen-Cooper-Schrieffer (BCS) reduced 
Hamiltonian 



^^BCS = E'^k("-kT +^-ki - 1) 
k 

- ^y^! 4'T'^-kM'^-k.l.CkT 
k,k' 



(122) 



where the (canonical) fermionic modes (electrons) are defined 
on a momentum space lattice with momentum index k and 
n\L(y ~ Cko-Cko- (with cr [). Again, this is a case where the 
dimension of the local Hilbert space is = 4 with a possi- 

and it turns out it is 



ble basis {|0),4^|0), 
exactly solvable. 

It is clear from this Hamiltonian that if the total charge of 
the pair (k "f , — k [) is one (i.e., the pair state is occupied by 
a single electron, rikt ~ 1 and n_ki = 0, or fiki = and 
n_ki = 1) it will remain equal to one after appUcation of 
i^BCS (Pauli blocking). We will show now that this conser- 
vation is derived from an SU{2) gauge symmetry of the BCS 
reduced Hamiltonian which becomes explicit when we write 
i/ecs in terms of the generators of su{2) ^ su{2). This con- 
stitutes another example of the use of alternative languages to 
unveil hidden symmetries and, to our knowledge, this has not 
been revealed in the way we will present below. 

We just need to note that the following sets of operators 















tl = t 




i3 = t/2 


















-2t 



1 



FIG. 9: Hopping processes for the bosonic Hubbard model written in 
the g-particles language, Eq. <12H . There are three different hopping 
processes ti, t2, and t^. Note, however, that the non- vanishing ma- 
trix elements of the kinetic energy Ta/X are the same (-2t) regardless 
of the hopping process. 



D. BCS reduced Hamiltonian 



'k 



(nkT + n-ki - 1) 



''kT'--ki 
-kiCkt , 



c 



(123) 



and 



5k 

•^k 
5,7 



*-kT -ki ' 

'^-ki'^kT ' 



(124) 



satisfy the spin su{2) commutation relations (fi, v,X = 

x,y, z) 



[5^, 5k/] 
K,rk^] 
?k" 



i<^kk'e^iiyA5k , 

■ s: A 
JOkk'e^iiyATk , 

. 



(125) 



The last set of local su{2) generators is associated to the 
conservation of the parity of the total charge in each pair 
(k t, — k [). Since these are local transformations, the sym- 
metry generated by them is an SU{2) gauge symmetry. To see 
this we just need to re-write i^ecs in the pseudospin language 



These general mappings between languages are not re- 
stricted to real space lattices, indeed one can find mappings 



^Bcs =2}^ekT^ 

k 



k,k' 



(126) 
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Since -ffBCS is only a function of the set of local su{2) genera- 
tors given by Eq. ( I123> it commutes with the set of generators 
given by Eq. ilM . 

In this way, this SU{2) gauge symmetry splits the to- 
tal Hilbert space into invariant subspaces labelled by the set 
of momenta A/bp which have single occupancy of electrons. 
Within each subspace (which corresponds to a depleted mo- 
mentum space lattice) the operators = c|^|clj^| are the 
single-flavor hard-core bosons whose algebra was defined in 
Eq. (|6}, and the Hamiltonian becomes 

k k,k' 

where k and k' mean sites in a momentum lattice k, k'^ Mhp- 
(This model is exactly solvable.) Thus, a superconducting 
BCS state is a Bose-Einstein condensate, in momentum space, 
of hard-core bosons (i.e., spinless JW particles). 



We have also seen in section IV^ that the su{N) genera- 
tors can be represented in terms of Schwinger-Wigner bosons. 
Therefore, we can use Schwinger-Wigner bosons to get a 
mean-field solution of the Hamiltonian (Eq. ( I129» . In this 
way we see how the change of language in the bosonic t-J 
model opens the possibility to get a simple and original solu- 
tion ileoli . 

These examples illustrate the process of unveiling hidden 
connections between different problems just by changing the 
language. We also saw in the bosonic Hubbard model exam- 
ple that the same process which unveils hidden symmetries re- 
duces the whole space to a direct sum of invariant subspaces. 
The effective number of degrees of freedom in each subspace 
is lower than the original one. Consequently, we can use a new 
class of languages ^ A F with a lower dim T to describe the 
original problem. This reduction not only simplifies the res- 
olution of the problem, but also establishes new connections 
with other physical systems. 



E. Bosonic t-J model 



In the previous section fVAl we introduced an su{N) lan- 
guage for the hard-core bosons (see Eq. ( I48» . In other words, 
we showed that the generators of su{N) in the fundamental 
representation can be used to describe a system of hard-core 
bosons with Nf = N — 1 different flavors. To illustrate 
a possible application of this result, we will consider here 
the case of spin 5'= 1/2 hard-core bosons (i.e., Nf = 2 and 
D ^ N = 3). In particular, we will show that the bosonic t-J 
model m 



(i.j> 



(128) 



with Sj = ^bj^Capb^p (cr denoting Pauli matrices), can be 
rewritten as an anisotropic SU (3) Heisenberg model. To this 
end we just need to use the identities given in Eq. ( I48t to re- 
place the bosonic operators by the corresponding su(3) gen- 
erators iS^^Q) to get 

Ht-j - E'^A^'-5'"'(^)'5'''(j) +^E'5""(j) ' (129) 

(iJ> J 

wifli ( J^^ = J^t,) 



Jni 



J 



Jw — J20 — t, 



J12 = Jii = J22 = i, B = ^l+^, (130) 
2 

and z is the coordination of the lattice. In this way we have 
established an exact mapping between the t-J model and an 
anisotropic SU (3) Heisenberg model with an applied 'mag- 
netic field' B. 



F. SO (3) Heisenberg Model 

As it was shown in Eq. j55L the generators of so(3) are 
the components of an antisymmetric tensor of order three. By 
contracting this tensor with the completely antisymmetric one 
Eq^i/, we can see that the three generators of so(3) transform 
like a vector L 



(131) 



It is well known that the so(3) and su{2) algebras are iso- 
morphic. As a consequence of this, we can establish a simple 
mapping between L and the three generators of su(2) S in the 
spin one representation 



S(j) = *L(j) 



(132) 



Using this basic result of group theory, we can map the 5*0(3) 
Heisenberg model 



(133) 



(iJ> 



= -UY^sfsf + J E (bl'bl ' blk. + H.c.) 

(iJ> (i.j)'^ 

into the S = 1 SU (2) Heisenberg model: 



S0(3) 
Heis 



SU(2) 
Heis 



-jES(i)-S(j). (134) 

(iJ> 



In other words, the S ^ 1 SU{2) Heisenberg model is equiv- 
alent to a model consisting of a gas of 5 = 1/2 hard-core 
bosons with a magnetic interaction along the y-axis and a su- 
perconducting term. 



VII. BROKEN VERSUS EMERGENT SYMMETRY 

Matter organizes in different manners depending upon the 
nature of its constituents and interactions. Symmetry and 
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topology are fundamental guiding principles behind this orga- 
nization. A symmetry transformation is a modification of the 
observer's point of view that does not change the outcome of 
an experiment performed on the same system. In mathemati- 
cal terms it is a transformation that takes the Hilbert space of 
states H into an equivalent one. Wigner's theorem asserts that 
any transformation T which conserves the transition probabil- 
ity between rays in H 

|(f^*|f$)p = (135) 

can be represented by a linear unitary or antilinear antiuni- 
tary map O {O^ = O^^) on H. Since the operation of time- 
reversal is one of the few relevant examples in physics which 
involves an antiunitary operator, we will only consider unitary 
mappings in the following. 

Symmetries may be classified as external or space-time 
(e.g., the Poincare group) and internal. The latter refers to 
the set of transformations that leaves the Hamiltonian of the 
system H invariant; i.e., these are the symmetries of the phys- 
ical laws. This set forms a group which is named the internal 
symmetry group Q and is defined as 

Q = {da} , with group elements ga, that satisfy [H, Qa] = , 

(136) 

with e representing the identity element, and where the num- 
ber of elements defines its order, which may be finite, denu- 
merable infinite (discrete), or non-denumerable infinite (con- 
tinuous). In general, groups of symmetries in physics are ei- 
ther finite or Lie groups (non-denumerable infinite). Besides, 
the group Q may be Abelian (i.e., [go,, go,'] = 0,Va,a') or 
non-Abelian, and local (also called gauge, meaning that the 
symmetry applies to subsystems of the original physical sys- 
tem) or global. Invariant physical observables, O, are those 
Hermitian operators which remain invariant under the sym- 
metry group Q, i.e., [0,gQ,] = 0. Every observable which is 
a function of the groups elements ga is a constant of motion 
since it commutes with H. Table HH shows representative ex- 
amples of physical models displaying different kinds of sym- 
metries. For each group element Qa there is a unitary operator 



Symmetry 


Discrete 


Continuous 


Global 


Ising (Z2) 
XYZ (D2h) 


classical XF(0(2)) 
Heisenberg (5'I7(2)) 


Local 


Gauge magnet {Z2) 


Gauge magnet ((7(1)) 
BCS (hidden SU{2)) 



TABLE IL Examples of models displaying different kinds of sym- 
metries. The group (or subgroup) of symmetries involved is written 
in parenthesis. 

(see above) that will be denoted as Oa = 0{Qa) which maps 
Ti, into an equivalent Hilbert space. The set {Oa} forms a 



representation of the group Q. A representation is an homo- 
morphic mapping of the group Q onto a set of linear operators 
O such that: 0(e) = 1, and OaOp = 0{QaQp)- The dimen- 
sion of the representation, dim(C'), is the dimension of the 
(vector) space on which it acts. By a representation we will 
mean a non-singular (in particular, unitary) dim(C') x dim(C') 
matrix representation. A representation is irreducible if its in- 
variant subspaces under the action of all the elements of the 
group are only and the full space. A completely reducible 
representation can be written as a direct sum of irreducible 
representations (irreps). The eigenstates of H that have the 
same eigenvalue £"„ form an invariant subspace 

HO^\^„) = 0„i/|*„) =^«Oa|*„) (137) 

meaning j'l'Jj) = Oa\^n) is also an eigenstate with the same 
eigenvalue. When the dimension of this invariant subspace 
is larger than one, the energy eigenvalue En is degenerate. 
The dimension of the degenerate subspace is equal to the di- 
mension of the representation of Q associated with the eigen- 
state l^'n). If the group Q is Abelian all the irreps are one- 
dimensional and there is no degeneracy induced by Q. 

Lie groups play a fundamental role in physics. There is 
a notion of continuity or closeness imposed on the elements 
of the group manifold Qa such that a finite transformation of 
the group can be generated by a series of infinitesimal ones. 
For a one-parameter continuous group the representations (the 
homomorphic mapping must be continuous) of its elements 
can be written 

0^{e) = exp[i6'X„] , (138) 

where 6* is a continuous parameter and X^'s are the generators 
of the Lie algebra. The representations of the group elements 
are defined such that ~ Q represents the identity operator II 
and an infinitesimal transformation 56 is expressed as 

(Da{5e) = t + i5e , (139) 

where the generators form a Lie algebra, 

[X^,Xp\=iCli3X^ , (140) 

with C'^p representing the structure constants of the algebra. 
Notice that the generators themselves are conserved quanti- 
ties; i.e., [H, Xa] = 0. 

Let us provide an example to show how these ideas are ap- 
plied. Suppose we have the following model Hamiltonian rep- 
resenting interacting spinless fermions 

(iJ> (ij> 

where (i, j) represents nearest neighbors in an otherwise bi- 
partite lattice (i.e., the union of two interpenetrating sublat- 
tices A and B). Among the elements of Q, there are two 
Abelian symmetries: one continuous and global U{\) related 
to charge conservation, and another discrete and local (stag- 
gered) Z2 related to a particle-hole transformation. The con- 
tinuous symmetry is realized by the unitary mappings 

=exp[*0E^j]' Ooc\ol = eMiO]c], (142) 
j 
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while the discrete one is realized by the identity and the 
particle-hole transformation 



Cj sublattice A 
-Cj sublattice B 



(143) 



where Jjs is one if j belongs to sublattice B and zero other- 
wise. 

In some instances the states of matter display the symme- 
tries compatible with the quantum equations of motion (i.e., 
symmetries of H) 



H\^[t))^idt\^[t)) , 



(144) 



while other more interesting situations are characterized by 
states with less (broken-symmetry scenario) or with more 
symmetries (emergent symmetry scenario). In subsection 
IVll BI we will introduce and expand on the latter concept. 



A. Broken Symmetry 

To put the concept of emergent symmetry in context let us 
start summarizing the, in principle independent but very pow- 
erful, broken-symmetry scenario. The broken-symmetry phe- 
nomenon, which is manifested by a lowest-energy state not 
having the full symmetry of the Hamiltonian Q but less, has 
been beautifully described in P. W. Anderson's Basic Notions 
of Condensed Matter Physics book jH. Here we will simply 
restate the main known results. 

To say that the ground state of H, |5'o), is invariant under 
Q means that 



Oa|*o) = |*o) , Va 



(145) 



If 1 4'o) is not invariant under a given symmetry operation Oa, 
we say that the symmetry is spontaneously broken. (Notice 
that if such a symmetry does not exist and jvPo) is not invari- 
ant that means that the symmetry is explicitly broken.) The 
broken-symmetry state is also called the ordered state. In gen- 
eral, a subset of Q, Qbs^ contains the transformations that do 
not leave I^Po) invariant, while a residual symmetry subgroup 
Qr remains. 

To make a quantitative distinction between the symmetric 
and the broken-symmetry phases we need to introduce the 
concept of order parameter, representing the supplementary 
variable needed to describe the lower symmetry state. By def- 
inition, the order parameter is a physical quantity which is 
zero for the symmetric phase and non-zero for the broken- 
symmetry state. However, this simple requirement is still too 
general since there are many different quantities that can sat- 
isfy this condition. The most natural choice is dictated by the 
symmetry which is spontaneously broken. The order parame- 
ter is then a physical quantity that transforms like a non-trivial 
representation of the symmetry group. The adequate choice of 
representation depends on the particular problem under con- 
sideration. For instance, out of a broken SU{2) symmetry, a 



system may have dipolar (usual magnetic ordering), nematic, 
or more sophisticated multipolar orderings. Each transforms 
like a different representation of the SU (2) group. One of 
the most important experimental challenges is the design of 
physical measurements that distinguish between the different 
representations of the broken-symmetry groups. 

Spontaneous symmetry breaking can be associated with 
collective phenomena, which is the relevant case for matter 
organization, or with a trivial phenomenon that occurs when 
part of the one-particle spectrum of H is degenerate. The lat- 
ter scenario is not of interest here because it does not lead to 
the important concept of ergodicity breaking. As an example, 
consider the ground state of an odd number of noninteract- 
ing spin- 1/2 fermions. The ground state is not SU{2) invari- 
ant (i.e, it breaks that symmetry) but as soon as we add one 
particle to the system the SU{2) symmetry is restored. On 
the other hand, when the phenomenon is collective, sponta- 
neous symmetry breaking is in general related to ergodicity 
breaking meaning that for given initial conditions the equa- 
tions of motion, although symmetric, cannot connect (because 
of kinematic or dynamical reasons) states that would other- 
wise restore those symmetries. The phenomenon of ergodicity 
breaking is not exclusive to systems with a macroscopically 
large number of degrees of freedom (thermodynamic limit), 
it may happen in systems with a finite number of degrees of 
freedom as well (e.g., a single particle in a double well with 
an infinite barrier in between). The first case is a result of the 
existence of a multitude of inequivalent representations of the 
observables of an infinite system. In other words, it is a result 
of the non-commutativity of the limits 



lim liimJj(n,F)^ lim lim ^Jj(n,F), 



(146) 



where ft is the volume of the system, and F is the general- 
ized external field that linearly couples to the order parameter 
-0. Notice that when the symmetry breaking is induced by 
ergodicity breaking, the subspaces which get disconnected in 
the dynamics of the system are still connected by symmetry 
operations of Gbs- This means that ergodicity breaking in a 
system with a Hamiltonian invariant under Qbs requires the 
ground state to be degenerate. 

Our last observation raises the following question: Is it 
possible to have spontaneous symmetry breaking when the 
ground state is non-degenerate? The answer to this question 
is yes. Symmetry breaking implies that the ground state is not 
invariant under the action of Qbs- Non-degeneracy means that 
the representation associated to this state is one-dimensional. 
Let's consider an example of this situation in which Q is the 
Abelian group of translations. The system consists of an A^^- 
site ring of noninteracting spinless fermions in the presence of 
a uniform magnetic field threading the ring whose total flux is 



= -t^(c]exp[i 



-i0]cj), (147) 



The ground state of this system is non-degenerate and has a 
non-zero total momentum (non-zero current); i.e., it does not 
belong to the trivial representation of the group of translations. 
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As usual, mathematical concepts have physical relevance 
whenever there are observable consequences. So, what are 
the consequences of symmetry breaking? They are: 

1. Generalized rigidity and long-range order 

With the appearance of a broken-symmetry state an order 
parameter tjj emerges which represents a measure of the de- 
gree of asymmetry in the broken-symmetry phase. The gen- 
eral problem of how to explicitly define an order parameter is 
deferred for a later section, here we will simply assume that 
we know i/' and it is defined as the expectation value of some 
space-dependent local observable {order field) ip{x.). Since 
the state is translationally invariant (in the thermodynamic 
limit) it happens that 

lim (V3(x)V3(x')) V'^ ■ (148) 

|x— x'l — »oo 

In other words, the broken symmetry state carries long-range 
correlations in the order field. Conversely, long-range correla- 
tions implies a broken-symmetry state, therefore, long-range 
correlations is a necessary and sufficient condition for the ex- 
istence of an ordered state. P.W. Anderson has called this phe- 
nomenon generalized rigidity. In his words iQ it is an ''emer- 
gent property not contained in the simple laws of physics, al- 
though it is a consequence of them." 

Given a physical Hamiltonian it is not straightforward, in 
general, to determine whether its ground state is invariant 
or not. However, under certain conditions (e.g., low space 
dimensionality) one can certainly establish that there is no 
broken-symmetry phase and, therefore, no long-range order 
Those conditions constitute the hypothesis of the Mermin- 
Wagner-Hohenberg theorem. The idea behind this theorem 
is that (thermal or quantum) fluctuations may destroy long- 
range order, so the conditions of the theorem are directly tied 
to the reasons that may cause strong fluctuations to the or- 
dered state. Continuous symmetries, low space dimensional- 
ity, short-range (constituent) interactions cause stronger fluc- 
tuations. Thus at finite temperature, T > 0, short-range spin 
models with a continuous symmetry in d < 2 do not show 
spontaneous ordering. At T = 0, the existence of a gap in the 
excitation spectrum (of a system with a continuous symmetry) 
precludes long-range order, but the presence of gapless exci- 
tations does not necessarily imply this order (Note that the 
Mermin-Wagner-Hohenberg and Goldstone (see below) the- 
orems are two complementary aspects of the same fact, and 
both follow in a unified way from a clever use of the Bogoli- 
ubov's inequality l6lll .) 

2. New massless particles: Nambu-Goldstone modes 

As mentioned above, continuous symmetries play a special 
role since degeneracies are non-countably infinite. If a con- 
tinuous symmetry is spontaneously broken, the spectrum of 
H generally has gapless collective excitations (soft modes). 
These emergent excitations (quasiparticles), which can be in- 
terpreted as new particles with zero mass, are called Nambu- 
Goldstone modes or bosons. They are the quantized excita- 
tions associated with a spatial twist of the order parameter 
The energy of the twist is typically proportional to the inverse 
of the system's length, i.e., it increases linearly with the wave 
vector k of the twist. (The energy of the resulting state van- 



ishes as k ^ 0.) The number of Nambu-Goldstone bosons 
is at most dim Q — dim Q^. The Anderson-Higgs mecha- 
nism in gauge theories provides an exception to this theorem, 
i.e., some (or all, depending upon the symmetry group of the 
gauge theory) Nambu-Goldstone bosons do not emerge even 
though the continuous symmetry is broken. The idea behind 
this mechanism is that after the symmetry is broken, the cou- 
pling of our ungauged system to a massless gauge field gener- 
ates a mass for the gauge field giving rise to massive bosonic 
excitations. 
3. Topological defects 

The breakdown of long-range order (with temperature, for 
instance) carries the formation of defect structures such as 
vortices and domain walls, each characterized by the type of 
singularity in the order parameter The topological stability of 
these defects is defined by the homotopy class of the manifold 
where the order parameter Uves, and work has been developed 
to mathematically classify these defect structures. The book 
by P.W. Anderson ||8i| provides an excellent introduction to the 
subject. 



B. Emergent Symmetry 

In the previous subsection, we revised the fundamental as- 
pects and the deep physical consequences of the concept of 
broken symmetry. We will introduce now another notion 
which complements the previous one and also plays a central 
role in the description of physical systems. This is the notion 
of emergent symmetry. In broken symmetry phenomenon, the 
symmetry of the considered system is lowered below a critical 
temperature. But is it possible to have the opposite situation 
where the symmetry of our system increases when the temper- 
ature is lowered? In a sense that will become clear below, the 
answer to this question is yes. However, we will see that the 
emergence of a new symmetry is not signaled by a phase tran- 
sition but by a crossover between two different energy scales. 
This new concept is useful to characterize the relevant degrees 
of freedom that dominate the low-energy physics and the na- 
ture of the quasiparticles that result from the interactions be- 
tween these degrees of freedom. 

The isolation of the relevant degrees of freedom for the de- 
scription of complex systems plays a central role in physics. 
In general, this process results from a careful separation of 
the different energy scales involved in the Hamiltonian H un- 
der consideration. The elimination of the irrelevant degrees of 
freedom can be achieved by different methods like perturba- 
tion theory, algebraic approaches, or renormalization group. 
The result is a new low-energy effective Hamiltonian, Heft, 
that is a valid description of the physical system below some 
characteristic energy Ec. The Hilbert space A^o of -ffgff is 
then a subspace of the Hilbert space H of the original Hamil- 
tonian H. Sometimes the spectrum of iJeff and the low-energy 
spectrum of H are exactly the same. However, in the more 
general case, the difference between both spectra is finite and 
can be made arbitrarily small by increasing the distance be- 
tween the corresponding energy scales in H. 

We will say that a new symmetry emerges whenever there 
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is a symmetry group G of transformations (which is not a 
group of symmetries of H) with elements of the form ~ 
X^k '^k n7=i ^^k is a c-number and rik an integer), where 
the unitary operator Uj acts on the local Hilbert space Tij , and 

[ga,Heff] = for any e G . (149) 

The transformations in G are defined in the original Hilbert 
space H and leave the subspace Mq invariant. [H = Tij 
is an arbitrary decomposition of H.] The condition = 
Xk'^knj=i means that the g^'s are linear combinations 
of nk-local operators (product of rik operators each of them 
acting on the local space Tij). In this way, we are excluding 
non-local transformations. This is because for each that 
commutes with H^ff, it is always possible to find a non-local 
operator that commutes with H and is identical to g^ when 
restricted to the subspace Aio- Note that if iJgff provides an 
exact description of the low-energy spectrum of H, we say 
that the emergent symmetry is exact. In contrast, if H^ff is 
only a very good approximation to the low-energy spectrum 
of H, the emergent symmetry is only an approximate concept. 
In other words, suppose that H^ff + H'^^f is the Hamiltonian 
that reproduces the exact low-energy spectrum of H. If, for 
instance, we have derived iJgff within perturbation theory to 
order n in the small parameter jy (ratio between the small and 
the large energy scales), H'^ff is of order 7y"+™, where m is a 
positive integer. Therefore, 

[5« , H,n+H'^f,] = [g„ , H^f,] = O (r;"+" ) for any .g„ € G . 

(150) 

Although in these cases the emergent symmetry is an approx- 
imate notion, the concept is still useful to identify the nature 
and properties of the low-energy quasiparticles. This is illus- 
trated in the last examples of this section. 

Let us start considering cases in which the notion of emer- 
gent symmetry is exact. During last few decades, exact dimer 
ground states were found for different quantum spin models 
1 10, 62]. The most famous example is the Majumdar-Ghosh 
solution of the Heisenberg spin- 1/2 chain with nearest and 
second-nearest neighbor interactions Ji and J2, respectively. 
The common characteristic of these quantum spin models is 
the emergence of a local U{1) gauge symmetry that gives rise 
to the conservation of each local dimer In general, the dimer- 
ized ground states can be written as 

1*0) = (g) IMh.h)) , (151) 

(ii42) 

where (ii, 12) denotes a pair of sites and each site belongs to 
one and only one pair. The wave function i2)) is the 

singlet state for the pair (ii,i2). Let us now introduce the 
following local U{1) transformation 

f(ii,i2) = expK2P(ii,i2)], (152) 

where P(ii,i2) = i2))(0s(ii, 12)! is a local projec- 

tor on the singlet state of the bond (ii,i2). It is clear that 
T(ii, i2)|^'o) = cxp[iai2]|^'o)- The generators of this U{1) 



gauge group are the projectors P(ii, 12) which count the num- 
ber of singlet states on each bond (ii, 12). Then, the set of 
bonds (ii , 12) provides the natural subsytem decomposition of 
H. It is important to note that in most of these cases the emer- 
gent symmetry is only present in the ground state which is sep- 
arated from the lowest energy excitations by an energy gap; 
i.e., Aio is a one-dimensional subspace. Nevertheless, this 
emergent property of the ground state determines the nature 
of the low-energy excitations, which are local spin triplets on 
each bond, that propagate on a sea of singlet dimers. For mod- 
els like the one of Eq. ( I88t or Ref. 1 13i], there is an invariant 
low-energy subspace given by a dimerized ground state and 
an arbitrary number of triplet excitations, with the condition 
that triplets cannot be created on dimers which are nearest- 
neighbors. The action of the Hamiltonian restricted to this 
subspace is invariant under a C/(l) gauge transformation and, 
again, the exact low-energy theory has an infinite number of 
symmetries not present in the original models. 

What is the origin of the emergent U{1) gauge symmetry 
of quantum dimer ground states of the form of Eq. dlSll l? 
For most of the cases, if not all, the origin is geometrical 
frustration. The particular connectivity of the lattice (or lat- 
tice topology) produces negative interference between the dif- 
ferent links that connect two dimers. When the interference 
is complete, the dimers become disconnected and the gauge 
symmetry emerges. This principle is by no means restricted 
to quantum dimer systems. We can imagine different local 
structures (instead of dimers) that get disconnected due to the 
same reason. For instance, currents localized in plaquettes 
can also result from an emergent SU(2) gauge symmetry in 
frustrated lattices. As shown in Ref. 16311 . this occurs for par- 
ticular fillings of the Hubbard model on a family of frustrated 
lattices. The guiding principle is always the same. In a certain 
region of parameters, local currents with two possible orien- 
tations become conserved quantities. The local chirality can 
be described with a pseudospin-1/2 variable that is the rele- 
vant degree of freedom to build a low-energy effective theory. 
Any weak physical interaction that breaks the emergent gauge 
symmetry can introduce a finite coupling between the local 
currents and produce different orderings of the chiral degrees 
of freedom. In the case of Ref. i63fl . the interaction that breaks 
the gauge symmetry is an inter-site Coulomb repulsion. This 
interaction induces an XF-like ordering of the local currents. 
The low-energy excitations are chiral-waves that are described 
by magnons in the pseudospin language. Following the same 
strategy, one can find physical Hamiltonians that give rise to 
different and, sometimes, unusual low-energy degrees of free- 
dom (see the example of Eq. ( I157t below). 

However frustration is not a requirement to have an exact 
emergent gauge symmetry. For instance, the Fermi liquid or 
the band insulator are two examples of emergent symmetries 
that do not involve any geometrical frustration. These phases 
of matter are not characterized by any broken symmetry since 
their corresponding ground states are non-degenerate. Ac- 
cording to Landau's theory of Fermi liquids |64|, the quasi- 
particles of the Fermi liquid have an infinite lifetime if they 
are right at the Fermi surface. This means that the low-energy 
spectrum of the Fermi liquid is invariant under local U{1) 
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transformations of the form cxp[i(/)kpcj^^C[j.^] (where cj^.^ is 
the creation operator of a quasiparticle with Fermi wave vec- 
tor lip)- The band insulator can also be characterized by a 
local U{1) emergent symmetry. In this case, the U{1) sym- 
metry can be factorized in real space and is generated by the 
local transformation: exp[iacjcj]. This is the mathematical 
expression for charge localization that characterizes the insu- 
lating state; the charge is locally conserved. Like in the case 
of the spin dimers, only the ground state of the band insula- 
tor exhibits the property of emergent symmetry. In contrast, 
the dimension of A^o for the Fermi liquid phase is equal to 
the number of wave vectors kp that are on the Fermi surface. 
Note that the Fermi surface is the manifold associated with the 
group of emergent symmetries whose topology characterizes 
the universality class of the Fermi liquid. 

Exactly and quasi-exactly solvable models provide also ex- 
amples for exact emergent symmetries. A model is quasi- 
exactly solvable when only part of the spectrum can, in a 
purely algebraic form, be exactly diagonalized. Let us call 
So the subspace generated by the exactly solvable part of the 
spectrum. Since H is exactly solvable when restricted to So, 
there is a set of operators that commute with H : So ^ Sq. 
If So is also the lowest energy subspace and the operators ga 
can be factorized as ga ~ J2k 11]=! ^^j' we have another 
case of emergent symmetry. As an example of a quasi-exactly 
solvable model that has also an emergent symmetry, we will 
consider the t-J^ chain |i6|]. The lowest energy subspace of 
this model can be mapped into the exactly solvable S=l/2 
XX Z chain. This means that H^_j : So ^ So has an in- 
finite number of symmetries 17,1 that are linear combinations 
of n-Iocal operators (products of n-body spin variables) l6^ . 
These are the quantum symmetries that make the XX Z chain 
an integrable problem. The relevant low-energy degrees of 
freedom of H^_j are holes which are attached to an anti- 
phase domain for the staggered magnetization. For < 8t, 
the system is a Luttinger liquid of particles that carry both 
electrical and topological charges. 

We will consider now the cases in which the emergent sym- 
metry is an approximate, albeit important, concept. In the sim- 
ple example that we describe below, a global SU{2) invari- 
ance emerges at low energies. The model is a d-dimensional 
hyper-cubic Kondo lattice with an attractive (U > 0) Hubbard 
interaction for the conduction band, an anisotropic Kondo in- 
teraction between the magnetic impurities and the conduc- 
tion electrons, and a Heisenberg antiferromagnetic interaction 
(J > 0) between the localized spins (cr =t, j) 

i i 

+ JifE(7^f + S[sf + Sfs]) + jE Si • Sj,(153) 

i (ij> 

where Si is the spin operator for the locaUzed moment at the 
site i, and = 1/2 X^^.r' 4r<r'Cir' with ly = {x,y,z\. 
The symmetry group of H^'^ is C/(l) x U{\). The corre- 
sponding generators or conserved quantities are the total num- 
ber of particles n\ and the z-component of the total spin 
X;(5'? + sf ). This model describes the competition between 



a Kondo-like system {Jk ^ t^) and an s-wave superconduc- 
tor coexisting with antiferromagnetic ordering of the local- 
ized magnetic moments {Jk ^ V'). We will consider here 
the large IJ limit U ^ Jk , t. In this limit, the low-energy 
subspace Aio is generated by states in which the sites of the 
conduction band are either empty or double occupied. In other 
words, the conduction states are non-magnetic and the Kondo 
interaction is therefore quenched. The low-energy effective 
model in A^o is 

= jE Si • + ^~ E (^^^j + ^l^i) + E 

- AE"! ' 
i 

where t — 2t^ /U, and jl — 2fi + zt. The hard-core bosons 
represent the local Cooper pairs 

The first observation is that the localized spins and the con- 
duction electrons are decoupled in H^f^. The original Kondo 
interaction is suppressed by the competing U term. As a con- 
sequence, the symmetry group of H^f^ is U{1) x SU (2). The 
U{1) symmetry is again associated with the conservation of 
the number of particles that now are Cooper pairs. The emer- 
gent spin rotational SU (2) invariance is explicit from the ex- 
pression H^f^ of Eq. jl54t . This is a simple example of a 
global SU{2) symmetry that only appears at low energies. 
The higher-order terms in Jk/U will remove this SU{2) 
invariance leaving the original U{1) symmetry of rotations 
around the z-axis. This means that one of the two Goldstone 
modes associated to the spontaneously broken SU{2) symme- 
try will acquire a small mass of order [Jk/UY with 71 > 2. 

In the same way an emergent global symmetry is helpful 
for identifying the nature of the quasiparticles, an emergent 
gauge symmetry provides a guiding principle for identifying 
the relevant degrees of freedom at low energies. In general, 
a noninteracting theory is characterized by local symmetries 
that express the independence of each particle. For instance, 
the translation of only one particle is a symmetry for a nonin- 
teracting gas. The inclusion of the interactions removes this 
local symmetry and the many-body problem becomes, in gen- 
eral, non-trivial. However, in many cases it is possible to find 
another locally gauge-invariant limit for the interacting prob- 
lem. In this case, the local gauge invariance emerges only at 
low energies and signals the appearance of effective degrees 
of freedom that become decoupled. For instance, the origin of 
antiferromagnetism is more transparent in the strong coupling 
limit of the half-filled Hubbard model than in the weak or in- 
termediate coupling regimes. Before giving a formal expres- 
sion for this statement, it is convenient to illustrate its meaning 
with this simple example. Let us consider a repulsive Hubbard 
Hamiltonian at half-filling (U > 0) 

ii^Hubb = E ^^a^^i'^ + CjV'^i'r) + ^ E ^'T"U 
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In the infinite- [/ limit, the low-energy subspace A^o is gen- 
erated by states having one particle per site, i.e., there is one 
spin 5 = 1/2 localized on each site. All states in the mani- 
fold Mq have the same energy, and this massive degeneracy 
is associated with an emergent SU{2) gauge symmetry. In 
other words, in this limit the local spins get decoupled and are 
free to rotate without changing the energy of the system. This 
tells us that the natural degrees of freedom to describe the sys- 
tem in the strong coupling limit are localized spins, instead of 
itinerant fermions. The presence of an emergent gauge sym- 
metry indicates that there is a limit in our original interacting 
theory for which some degrees of freedom become noninter- 
acting. The internal structure of these degrees of freedom is 
determined by the gauge symmetry group. If the system in 
our example is close to the gauge invariant limit, U ^ t, but 
finite, the relevant degrees of freedom are still the same but 
they become weakly interacting. This is the origin of antifer- 
romagnetism in the strongly interacting Mott insulators. Since 
the large-C/ limit of the Hubbard model is very well known, 
it may seem to the reader that the use of the new concept of 
emergent symmetry just provides a complicated way to de- 
scribe a simple phenomenon. However, this is not the case for 
the non-trivial examples that are described below. In addition, 
we will see that the concept of emergent gauge symmetry is 
important as a guiding principle to find new states of matter 
that result from strongly interacting systems. One has to keep 
in mind that even antiferromagnetism remained as a hidden 
phase for a long period of time. 

Let us consider the case of a gauge symmetry group that is a 
direct product of local symmetry groups, G ~ (^j^ Gk, where 
each local group Gk acts on the local space Tik- If the reduced 
Hilbert space A4o admits the decomposition A4o = {^j^ Hk, 
this will provide the natural basis to write down Hgff. In par- 
ticular, if Gk is a continuous group, the generators of Gk 
are conserved quantities at low energies because they com- 
mute with Hsff. These generators are physical degrees of free- 
dom and their conservation imply that they are noninteracting. 
Note that these effective degrees of freedom are the bricks to 
build new phases out of interacting systems. As soon as we 
move away (but not too far) from the gauge invariant limit, 
these degrees of freedom will interact, producing, in some 
cases, novel types of orderings. The main goal in the rest of 
this section will be illustrating this phenomenon with different 
examples. 

We will consider now the spin- 1/2 ladder of Fig. ^| de- 
scribed by the following Hamiltonian 

H^^ = jJ2 • Si.' +J' J2 ■ Si+1., (157) 

ijy.i/' i,z^— 1,4 

with J, J' > and I < < A. In the limit J' <C J, 

the low-energy subspace of H^^ only contains states in which 
each square plaquette i is in a singlet state. This is clear when 
we analyze the spectrum of an isolated plaquette. The energy 
of the eigenstates only depends on the total spin St- E{St) = 
St{St + i)J- Thus, the two possible singlets {St ~ 0) are 
the lowest energy states, with the S't = 1 and St = "2 states 
having eigenvalues 2 J and 6 J, respectively. 

We will classify the two local singlet states according to the 




FIG. 10: Four-leg spin ladder. In each plaquette i, there are 4 spin- 
1/2 which interact through a Heisenberg exchange term of strength J. 
The interaction between plaque ttes is J', and the links are indicated 
in the figure. 

eigenvalues of the reflection symmetry plane along the diag- 
onal of the plaquette (see Fig. One of the singlet states 
\Si,s) is symmetric under this reflection while the other one, 
jS'i.a), is antisymmetric. These two singlet states can be de- 
scribed with an effective pseudospin t = 1/2 variable. We 
will represent the symmetric state |S'i,s) = |t) by the eigen- 
vector of with eigenvalue 1/2 and the antisymmetric one, 
I'S'i.o) = ID- by the other eigenvector of r^. In the infinite- 
J limit, the low-energy degrees of freedom Ti are completely 
decoupled. This situation is similar to the infinite-C/ limit of 
the Hubbard model. However, in the present case, the pseu- 
dospin variable ri does not represent a magnetic degree of 
freedom. By means of a canonical transformation to second 
order in J'/ J, we can derive the following low-energy effec- 
tive Hamiltonian 

where Jeff = -{J'f/J- H^k is an Ising-like anisotropic 
Heisenberg model. Since Jeff < 0, the ground state is a 
fully polarized ferromagnet that has two possible orientations 
{Z2 broken symmetry). In terms of the original variables, 
the fully polarized states correspond to the direct product of 
the symmetric and the antisymmetric singlets: l^i.,,) and 
\Si.a)- The Z2 X U{1) symmetry of H^lf is an emergent 
symmetry because H is not invariant under the corresponding 
transformations. These symmetries are explicitely removed 
when higher order terms in J'/ J are added to H^^- These 
terms will stabilize one of the fully polarized ferromagnetic 
solutions. However, the presence of these terms does not mod- 
ify the nature of the quasiparticles of i?|ff . These novel quasi- 
particles, which we will name kekulons l66ll . are magnons in 
the pseudospin language whose spectrum is gapped due to the 
Ising-like character of H'^^. In the original language, these 
magnons are massive singlet waves. The main effect of the 
higher-order corrections to H'^^ is a renormalization of the 
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mass (gap) of these quasiparticles. Therefore, even though the 
emergent symmetry is approximate in the present problem, its 
most important physical consequences remain valid. 

We have seen that there is a systematic procedure for ob- 
taining new types of orderings. The procedure starts from 
some limit in parameter space, in general the strong coupling 
limit, in which effective low-energy degrees of freedom are 
simultaneously stabilized and decoupled. The interactions be- 
tween these degrees of freedom appear when we move away 
from this limit and the symmetries or approximate symmetries 
of these effective interactions have important consequences 
for the properties of the ground state and the low-energy exci- 
tations. These degrees of freedom are the bricks for the novel 
orderings that result from the effective interactions and sim- 
plify the description of the new phase. These novel orderings, 
that are in general associated with the presence of competing 
interactions, pose a challenge for experimental physicists to 
develop probes that are sensitive to the new phases. 



VIII. BROKEN SYMMETRY PHASE TRANSITIONS 

It is widely recognized that the resolution of interacting sys- 
tems at a microscopic level requires sophisticated techniques 
which go beyond the traditional approaches. Much of the cur- 
rent understanding on quantum phases of matter and their cor- 
responding phase transitions is due to renormalization group 
analysis or the existence of a few exact solutions. The lat- 
ter, the most desirable scenario, is also the hardest since one 
seldom encounters solutions in spatial dimensions larger than 
one. For strongly interacting systems the situation worsens 
because it is difficult, if not impossible, to identify a small pa- 
rameter: Their ground states are susceptible to different quan- 
tum orderings and low energy excitations because of com- 
peting interactions that abound in nonlinearities. Often these 
systems are near quantum criticality which makes their study 
even more challenging. 

A fundamental notion of universality (or equivalence) nat- 
urally emerges from our dictionaries. This notion refers to the 
fact that many apparently different problems in nature have 
the same underlying algebraic structure and, therefore, the 
same physical behavior In this way, as a result of unveil- 
ing the symmetry structure of the problem, there is a concept 
of physical equivalence hidden in those dictionaries. If the 
complete Hilbert space of the problem maps onto another in 
a different language, the concept of universality applies to all 
length and time scales. On the other hand, if only certain in- 
variant subspaces of the original Hamiltonian map onto an- 
other physical system, then, universality will only manifest at 
certain energy scales. 

It was Landau 0] who first noticed the crucial role of sym- 
metry for phase transitions. Since the symmetry of a given 
state cannot be changed continuously, different symmetries 
must be associated to states of matter which are qualitatively 
different. To characterize the symmetry of a given state. Lan- 
dau introduced the notion of an order parameter This param- 
eter is zero in the symmetric state and has a nonzero aver- 
age when the symmetry is broken. Another property of this 



parameter is that it is non-invariant under at least one of the 
transformations of the system Hamiltonian. As Landau rec- 
ognized long time ago, the order parameter is the relevant 
physical quantity to build a macroscopic theory of thermody- 
namical phase transitions. However, the search for the order 
parameter characterizing a broken symmetry state can be a 
highly non-trivial task. 

There are examples in nature of systems entailing hidden 
microscopic order parameters which do not correspond to 
macroscopic variables, i.e., they cannot be coupled to any ex- 
ternal physical field. To find out these hidden order parame- 
ters we cannot avoid a microscopic description of our physical 
system. It is important to notice the practical consequences of 
this search; A new microscopic order parameter amounts to 
predicting a new possible state of matter The glassy mate- 
rials are good examples of systems belonging to this class. 
Since the order parameters characterizing a broken continu- 
ous symmetry are associated in general to the generators of 
the symmetry group, one can imagine that a language based 
on symmetry generators can shed some light into the search 
of the possible order parameters for a given system. It is the 
purpose of this section to illustrate this idea with different ex- 
amples and to introduce the notion of hierarchical language. 
One of our goals is to show the fundamental role played by the 
hierarchical languages in the classification of orders in matter 



A. Classification of Order Parameters 

To illustrate the general procedure of getting and classi- 
fying the possible order parameters of a given model, we 
start discussing the family of Heisenberg models with SU{N) 
spins in the fundamental representation 

H^T^ ^jY,S'^-'{i)S'''^{i). (159) 

(A summary of the methodology to obtain and classify or- 
der parameters is presented at the end of section lVlIIBI ) For 
J < 0, the ground states of this family of Hamiltonians can 
be exactly computed and it is a fully polarized SU (N) ferro- 
magnet. The order parameter is the SU [N) magnetization 
'5'"' = which has iV^ — 1 (number of genera- 

tors of su{N)) components. If J > 0, and the ground state 
has antiferromagnetic long range order, the order parameter is 
the staggered SU{N) magnetization S^j = J2} e'^'^S'"'{i) 
(Q is the antiferromagnetic wave vector). By adding terms 
which break the SU [N) symmetry explicitly, we can get 
lower symmetry order parameters associated to the subgroups 
of SU{N). 

We can ask now what is the expression of these order pa- 
rameters in a spin SU (2) or in a bosonic representation of 

^Heis^^- The answer to this question will allow us to explore 
new phases which can appear in the spin or the bosonic equiv- 
alent models. In addition, we will find that the high symmetry 
points of these models correspond to situations where two or 
more different phases are coexisting. 
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We start with the spin SU{2) representations of -H^Heis^'- 
The local Hilbert space of this SU (N) Heisenberg Hamilto- 
nian has dimension N so we can use an SU{2) spin S = 
{N — l)/2 to represent this model. In this representation the 
Hamiltonian Eq. ( I159l l has the following form 



SU(N) 
Heis 



JV-1 



(ij> '=0 



(160) 



with the values of a/'s determined from the following system 
of equations 



ai{[l ~ S{S + 1)]' - + 1)]'} - 2(-l)^ , 

Y.^A\yx - s{s + 1)]' - + 1)]'} = , 

=0 

"'{[^2 - S{S + 1)]' - [1 - S{S + 1)]'} = , (161) 



!=0 



1=0 



where = m{2m + 1), 3^2 = n{2n — 1) and m, n are 
integers. When N is an odd integer (i.e., S is integer) m and 
n satisfy < ?7i < 5, and 1 < n < S", while for even N 
(i.e., S is half-odd integer), < ?7i < 5* - 1/2, and 1 < rt < 
S + 1/2. It is easy to check that there are — 1 linearly 
independent equations and N variables ai . The undetermined 
variable corresponds to an additive constant. 

The simplest case corresponds to SU{2). Here we get the 
spin 1/2 (S'=(2-l)/2) version 



i^^f)=2j5:s..sj 



(162) 



(iJ> 



To understand the general procedure we continue with the first 
non-trivial example, i.e.jSU{3). This particular example has 
been analyzed in Ref. 14611 . By solving the above system of 
equations we get that the spin one (S'=(3-l)/2) version of the 
SU{3) Heisenberg Hamiltonian is (up to a constant) 



SU(3) 
Heis 



jV[Si.Sj + (Si.Sj)2] 



(163) 



To fix ideas let us start with the simple ferromagnetic case 
J < 0. As mentioned above, the ground state of the fer- 
romagnetic SU{3) Heisenberg model is the state with max- 
imum total SU (3) spin S. The order parameter associated 
to this broken symmetry is the total SU{3) magnetization 
gfj.i' _ S^'^d) which has eight independent components. 
We now raise the following question: what is the order param- 
eter of the equivalent model Eq. (I163> written in the {SU{2)) 
5=1 language? To answer this question we need to write down 
the relation between the components of 5^'' and the 5=1 gen- 
erators. The fundamental theorem guarantees the existence of 
these mappings. From Eq. (15 0> and the generalized JW map- 



pings (section fV C> . we get (the site index j is omitted) 



5^ 

sy 

{5^,5^} 
{5", 5"} 

{S\S'} 



1 



V2 



-1 



V2i 
v2 

_i(501_520^02_^10)^ 

V2i 



-5" 



(164) 



The first three operators are the components of the 5=1 ferro- 
magnetic order parameter, while the second five are the com- 
ponents of the spin-nematic order parameter (components of 
the bilinear symmetric traceless tensor). The traceless condi- 
tion implies that {Sv)^ - | = -((5^)2 - |) - ((5^)2 - |). 
In this way, we see that by rotating the SU (3) ferromagnetic 
order parameter it is possible to evolve from a ferromagnetic 
to a spin-nematic ground state and vice versa. This means that 
the 5t/(3) invariance of the spin Hamiltonian gives rise to co- 
existence of both phases. Since the algebra su{3) provides a 
fundamental representation when the local Hilbert space has 
D — 3, any local and linear operator can be written as a lin- 
ear combination of the identity and the su(3) generators. The 
5[/(3) magnetization is the highest-dimensional order param- 
eter for the 5=1 problem. By reducing the symmetry of the 
Hamiltonian to any subgroup G of SU (3) we can obtain lower 
symmetry order parameters. The general procedure consists 
in building a new basis for the generators where each opera- 
tor transforms according to an irreducible representation of G. 
For instance, if we take G = SU{2) we can do this classifica- 
tion by inverting the Eqs. ( I164> 



500 

^01 
^20 

5°2 
512 



5^(5^ + 1) 



1 

3 ' 



1 



2V2 
1 

2^/2 
1 

2V2 
1 

2V2 

i 



[s- 



[s- 



{5+,5-'}] 

{s-,s^]] 



[5+ - {5+,5"}] 



-{s\sy} + {s'-f + -{s'f^ 

i{5^5n + (5^)2 + 1(5^)2 
2i 2 



1 



(165) 



In this way we see that the ferromagnetic and the spin-nematic 
order parameters form a particular basis of su(3) generators 
in the fundamental representation. In other words, we can 
write down any linear operator as a linear combination of the 
components of both order parameters. Therefore, if a 5=1 



33 



LA-UR-OO-XXXX 



February 6, 2008 



Advances in Physics 



Hamiltonian is SU (2) invariant the local order parameter can 
be either ferromagnetic or spin-nematic since it can be written 
as a linear combination of both. This exhausts all possible 
order parameters for an SU{2) invariant 5=1 Hamiltonian. 

Following the previous procedure we can now classify the 
possible local order parameters of an SU{2) spin S Hamilto- 
nian. The ordinary magnetization is the only phase that can 
be derived from S=\I2 spins. It is well-known that a spin- 
nematic phase like the one described above cannot exist for 
5= 1/2. For S'=l, we have seen that there is an additional spin- 
nematic phase which appears in a natural way from the con- 
nection with the SU{?>) group. In general, we can write down 
a spin S Hamiltonian in terms of the generators of SU (25*+ 1) 
in the fundamental representation. To determine the possible 
local order parameters of an SU [N] invariant spin S Hamil- 
tonian we have to reduce this space of generators according to 
the irreducible representations of SU{2). It is easy to check 
that those representations correspond to the totally symmet- 
ric tensors of rank I < 2S, i.e., the possible order parameters 
are generated by application of the SU (2) transformations to 
the set: {5^ (5^)^, • • • , {S"")^^}. (S^'Y is the highest weight 
operator of the SU{2) representation associated to the totally 
symmetric tensor of rank I. The dimension of this representa- 
tion is 21 + 1. Therefore, the dimension of the space spanned 
by these operators is X]?fi(2^ + 1) = (25 + 1)^-1 which 
coincides with the dimension of the space of generators of 
SU{2S+ 1). In this way we see that there are 25 independent 
local order parameters for a spin 5 problem. The first two cor- 
respond to the local magnetization and the local spin-nematic 
order parameters. To our knowledge, there is no special name 
for the other multipolar orderings. 

This last example illustrates the general procedure to fol- 
low in order to determine and classify the possible local order 
parameters of a given system. It is important to remark that 
this classification can be made in any representation j43l . For 
instance, we can re-write the SU (3) Heisenberg Hamiltonian 
in terms of 5=1/2 hard-core bosons. In that language the lo- 
cal order parameters may be the local magnetization and the 
Bose-Einstein condensation order parameter In this way 
our algebraic procedure allows one to identify the possible 
local order parameters of a given system. In addition, the co- 
existence of more than one order parameter (more than one 
phase) can be described in a unified way by using an adequate 
language. 

Another connection, which is useful to illustrate the pur- 
pose of this section, is the one relating the SU (N) Heisenberg 
Hamiltonians to t- J-like Hamiltonians for hard-core bosons 
with spin 5 = (iV - 2) /2 (i.e., A^/ = - 1 different flavors). 
To this end we have to use the bosonic expressions for the 
generators of SU {N) (Eq. ( I48» introduced in section IVl Us- 
ing these expressions we can re-write the SU (N) Heisenberg 
Hamiltonian in the following way 

^SU(N) _ ^SU(N-l) , jV- r^t I +5t 5 1 

rV^ - - - 1 _ 

+ J 2^ - z J ^ 2^ : (166) 

(i.j> j 



where nj YULx "-j"' B^}"^ ^' is given by 

il^^}s~''^ = J E S^^iS)S''^{S). (167) 

This mapping is valid for any spatial dimension d. If we write 
down the same Hamiltonian using fermions instead of hard- 
core bosons, a gauge field appears in d = 2 due to the pres- 
ence of the operators which transmute the statistics. When 
d — \ the Hamiltonian is exactly the same for fermions and 
hard-core bosons (and anyons, in general). In this particular 
language we can check that the ferromagnetic SU (N) order 
parameter describes the coexistence of all the magnetic phases 
associated to the SU{N — 1) Heisenberg Hamiltonian (for in- 
stance, ferromagnetic and spin-nematic phases for A^ = 4), 
a Bose-Einstein condensate for each of the A^ — 1 different 
flavors, and a homogeneous (k = 0) charge density wave. 

Let us start analyzing the simplest N = 2 case (Eq. (I162t ). 
This case corresponds to spinless bosons with a kinetic en- 
ergy term and a nearest-neighbor density-density interaction. 
There are no spin degrees of freedom for this particular case 
and the SU (2) invariance is then associated to the charge de- 
grees of freedom. The two phases which coexist in this case 
are the Bose-Einstein condensate and the uniform charge den- 
sity wave order For A^ = 3 we have the usual t-J Hamil- 
tonian for spin 5=1/2 hard-core bosons (or fermions in Id). 
In this case the SU (3) symmetry gives rise to a coexistence 
between 5=1/2 ferromagnetism, Bose-Einstein condensation 
in both flavors, and the uniform charge density wave. These 
orders are different components of the same SU (3) order pa- 
rameter (see Eq. (I50» . In this particular language, = 5^°, 

foj = 5^°, B| = 5"^, and = 5"^ are the components of the 
order parameter for the spin up and down Bose-Einstein con- 
densates. For the magnetization we get = ^('^^^ ~ 5^^), 

that JT-j = § — 5°" can take any value by making SU (3) rota- 
tions of the ground states is another manifestation of the Bose- 
Einstein condensation. In addition, nj is the order parameter 
associated to the uniform charge density wave (if the SU (3) 
magnetization is oriented along the charge density wave di- 
rection we obtain nj = 1 in each site). For A^ = 4 we get 
a model which describes a gas of 5=1 hard-core bosons with 
the bilinear-biquadratic magnetic interaction of Eq. (( I163» . In 
this case we get coexistence of ferromagnetism, spin-nematic, 
Bose-Einstein condensation for the three different flavors, and 
a uniform charge density wave. As in the previous case there 
is an additional phase when A^ is increased by one. 

Therefore, this family of SU{N) Heisenberg Hamiltonians 
naturally describes a multi-phase behavior when the original 
language, based upon the generators of SU{N), is translated 
into another based on generators of a lower symmetry group. 
This illustrates a general rule which applies to any pair of 
groups (G, Q) such that G is a subgroup of Q, and the order 
parameter P transforms according to an irreducible represen- 
tation of Q. In general, this representation will be reducible 
under the operations of the lower symmetry group G, i.e., we 
will be able to express the given representation as a direct sum 
of representations which are irreducible under the application 
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of G. We can associate an order parameter to each of these 
representations. If we add the dimensions of each of these 
order parameters, we will get the dimension of P. Therefore, 
by using the lower symmetry language, we obtain different 
phases, and each of them is characterized by one order pa- 
rameter p-y . 

In the previous analysis we have not discussed the possi- 
ble spatial dependence of the order parameter In other words, 
we assumed that it was homogeneous over the entire lattice. 
However, there are many instances in nature where the or- 
der parameter is non-uniform. The antiferromagnetism is one 
of the most common examples. The previous analysis can 
be easily extended to the case of non-uniform order parame- 
ters with a well-defined wavevector k. We just need to rec- 
ognize that the non-uniform order parameters can be written 
as a Fourier transform of the local order parameter p(j) if 
the model considered has lattice translation invariance. In- 
deed, the presence of a non-uniform order parameter indicates 
that the translation symmetry has been spontaneously broken. 
Therefore the expression for the non-uniform order parameter 
p(k) is 



P(k) = ^^exp[ik-j] p(j) 



(168) 



Since the previous analysis is applied to the local order pa- 
rameters p(j), Eq. ( I168> shows that conclusions are trivially 
extended to non-uniform global order parameters. 



B. Hierarchical Languages: Tlie Quantum Phase Diagram of 
the Bilinear-Biquadratic Heisenberg Model 

We have seen that the local order parameter acquires its 
simplest form when it is expressed in terms of the hierarchical 
language. In addition, the generators of this language exhaust 
all possible local order parameters which may result from the 
solution of the problem under consideration. In other words, 
any local order parameter can be written as a linear combina- 
tion of generators of the hierarchical language. The Hamil- 
tonians considered above are special cases since they have an 
SU{N) invariance and therefore correspond to high symme- 
try points of an eventual phase diagram. We just considered 
those cases as the simplest examples of coexistence of differ- 
ent phases and unification of order parameters. We will show 
now that high symmetry is not a requirement for the success- 
ful appUcation of the present formalism to the determination 
of quantum phase diagrams. To this end we will consider 
the most general isotropic SU{2) S=l model with nearest- 
neighbor interactions in a hypercubic lattice (an overal J > 
factor is omitted in the following) 

= V2^[cos(/)Si • Sj +sin(/)(Si • Sjf] , (169) 

(iJ> 

a model already introduced in Eq. ( I83t . The parameter (f> 
sets the relative strength between the bilinear and biquadratic 
terms. Eq. ( I163> corresponds to the particular case = 57r/4. 
Indeed, as it is shown below, there are four isolated values of 



<j> = {7t/4, 7r/2, 57r/4, 3tt/2} for which is SU{3) invari- 
ant. The only symmetry which is present for any value of (p 
is the global SU{2) invariance since is a function of the 
scalar products Si • Sj. 

The Hamiltonian i/^ has been the subject of several studies 
in the last two decades [14, 67, 68, 69, 70, 71, 72, 73, 7i 
ItMIt^IttIIt^I . nevertheless, the complete characterization of 
the different phases was not, until now, completely solved. A 
semiclassical treatment for d > 1 fl4ll indicates that there are 
four different phases: the usual ferromagnetic (7r/2 < cj) < 
57r/4) and antiferromagnetic {3tt/2 < cj) < 7r/4) phases are 
separated on both sides by collinear- (57r/4 < cj) < 3tt/2) 
and orthogonal-nematic [tt/A < (f> < tt/2) orderings. We will 
show below that the collinear- and orthogonal-nematic phases 
obtained with the semiclassical approximation are replaced by 
uniform- and staggered-nematic orderings, respectively. 

As we have seen in the previous subsection, the SU{3) 
spins in the fundamental representation and the 5=1 SU{2) 
spins are two equivalent languages. In addition, we have 
shown in the sections FV Al and lV CI that the SU{3) spins and 
the 5=1 SU (2) spins can be respectively mapped onto 5=1/2 
hard-core bosons. We will use now these transformations to 
map the spin one Hamiltonian i?^ onto its SU{3) spin ver- 
sion. For pedagogical reasons, it is convenient to use the 
5=1/2 hard-core bosons as an intermediate language. 

In section IV^ we introduced a spin-particle transforma- 
tion connecting SU{N) spins and multiflavored hard-core 
bosons (JW particles, in general). In particular, the funda- 
mental (quark) representations of SU (N) were mapped onto 
an algebra of hard-core bosons with Nf ^ N — 1 flavors (see 
Eq. (I48» . For = 3 the hard-core bosons have two flavors 
(a =1, 1) which can be associated to an internal spin 5=1/2 
degree of freedom. A compact way of writing the SU (3) spin 
in terms of hard-core bosons is 



' 3 



b] 



b 



'^ji " 3 



(170) 



It is straightforward to write down each generator of the su(3) 
algebra in terms of the Gell-Mann (traceless Hermitian) ma- 
trices 



Ai 



A4 = 
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with the result 





Ai 


+ iA2 

' 
Z 


Hi - 


Ar 


— iAy 




2 


hi = 


A4 


+ tAs 




2 


1 


A3 


+ 'r 


" 3 


2 





Ai - zA2 



Ag + iAv 



2 

A4 



lAs 



A^ 
V3 



This constitutes the Cartan-Weyl representation of su(3) and 
illustrates the fact that a 5=1 SU{2) spin can be equally rep- 
resented in terms of quark fields. 

In the same way we wrote in Eq. jl70> the generators of 
SU{2i) in the fundamental representation, we can write down 
the corresponding expressions for the generators in the conju- 
gate representation 



50) = 



' 3 



-Hi 



(171) 



When the S=l operators are replaced by the corresponding 
functions of SU{?>) generators in the fundamental and the 
conjugate representations, it turns out that i/^, up to an ir- 
relevant constant, is a linear combination of the ferromagnetic 
and the antiferromagnetic 5[/(3) Heisenberg models 

= A/2^[cos05^''(i)5'^'^(j) 

(iJ> 

+ (sine/)- COS . (172) 

Repeated greek superindices are summed and the site index i 
runs over one of the two sublattices. This expression for iJ^ 
illustrates the very important result that any nonlinear inter- 
action in the original representation is simply a bilinear term 
in the new representation when mapped onto the highest rank 
algebra 14^ . In particular, as mentioned above, there are cer- 
tain special points in parameter space where the Hamiltonian 
is highly symmetric. For example, for = f and is 
explicitly invariant under uniform S'C/(3) transformations on 
the spins |'4S], while for = ^, is explicitly invari- 

ant under staggered conjugate rotations of the two sublattices. 
These symmetries are hard to identify in the original spin rep- 
resentation but are manifest in the SU (3) representation. 

In the following we will concentrate on the determination of 
the quantum phase diagram of H^p for spatial dimensions d > 
1 to avoid the strong effects of quantum fluctuations which can 
restore the continuous symmetry when d = 1. In the previous 
section we have analyzed the high symmetry point (j> = 
We found that the ground state has a non-zero order parameter 



(173) 



associated to a broken continuous SU (3) symmetry. This or- 
der parameter is the uniform SU (3) magnetization and corre- 
sponds to the coexistence of a ferromagnetic and a uniform 



spin-nematic ordering (see Eq. ( I164t ). This indicates that 
(j) — is a quantum phase transition point separating a fer- 
romagnetic phase from a uniform spin-nematic one. Let us 
consider now the related point (f> = j which differs in an 
overall sign from the previous case. This sign changes the in- 
teraction from ferromagnetic to antiferromagnetic. Therefore, 
for this new high symmetry point we expect to get a ground 
state characterized by the staggered order parameter 



(174) 



where Q is the antiferromagnetic wave vector It is clear from 
Eq. ( I164l i. that this staggered SU (3) order parameter corre- 
sponds to the coexistence of the staggered SU{2) magnetiza- 
tion 



MsT = E^^P[^Q-j] 



and the staggered nematic order parameter 



(175) 



(176) 



Nj is the symmetric and traceless component of the tensor 
obtained from the tensorial product of two vectors Sj. Hence, 
(j> ~ J is a transition point separating the usual antiferromag- 
netic ordering from a staggered spin-nematic phase character- 
ized by the order parameter of Eq. ( I176> . 

We will consider now the other two high-symmetry points, 
(j) = f , -Y"- For (p = the SU (3) symmetry is generated by 
the staggered operator 



j6S 



(177) 



where A and B denote the two different sublattices of a hy- 
percubic lattice. In this case, we have a ferromagnetic interac- 
tion between S{i) and S{j), and then S+ is the order parame- 
ter characterizing the broken SU (3) symmetry of the ground 
state. It is interesting to note that when the SU (3) order pa- 
rameter S+ is reduced with respect to the SU (2) group, the 
two coexisting order parameters are the staggered magnetiza- 
tion (see Eq. ( I175» and the uniform nematic order parameter 



N 



En. 



(178) 



In other words, if we apply an SU{3) rotation generated by 
S+ to the staggered magnetization we get the uniform nematic 
order parameter, and vice versa, the uniform nematic order pa- 
rameter is rotated into the staggered magnetization. This can 
be immediately seen by writing down the components of the 
local SU (2) magnetization and the nematic order parameter 
as a function of the local generators of su(3) in the conjugate 
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Comparing these expressions to the ones in Eq. ( I164l i. we see 
that when we change from to there is a change in 
sign for the three components associated to the magnetiza- 
tion, while the five components corresponding to the nematic 
parameter remain the same. Then, it is clear that 5+ describes 
the coexistence of a staggered magnetization and a uniform 
nematic ordering. Therefore, the conclusion is that </> = -^j^ 
separates an ordinary antiferromagnetic phase from the uni- 
form nematic ordering. 

The last high symmetry point to be considered is = f . 
In this case the coupling between and turns out to 
be positive, i.e., antiferromagnetic, and therefore we expect to 
get a broken continuous symmetry characterized by the order 
parameter 

S- =Y.S{S) -Y^SiS) . (180) 

From the considerations above, it is clear that 5_ describes the 
coexistence of ferromagnetism (uniform magnetization) and 
staggered nematic order Hence, </> = |- is a transition point 
separating these two phases. 
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TABLE III: Order parameters describing the different phases of the 
bilinear-biquadratic S=\ Heisenberg model for d > 1. indicates 
the phase boundary where the two phases in parentheses coexist. 

In this way, by identifying the high-symmetry points of B.^ 
we have determined the quantum phase diagram of this model 
(see Fig. II 1>. In addition to the transition points, we have ob- 
tained explicit expressions for the order parameters associated 



to each phase for any d > 1; these are summarized in table 
nm We can also predict from this analysis that the four tran- 
sition points (high-symmetry points) correspond to first order 
quantum phase transitions. In each phase, the corresponding 
order parameter has a a finite value and they coexist pairwise 
at the high-symmetry points. However, as soon as we depart 
from this point in one or the other direction in 0, the SU (3) 
symmetry is removed and one of the order parameters goes 
discontinuously to zero. In other words, the states with pure 
magnetic (ferro or antiferromagnetic) and nematic orderings 
belong to different representations of SU{2) (the remaining 
symmetry) so only one of them remains as the ground state 
when the SU{3) symmetry is lifted. 




FIG. 11: Quantum phase diagram of the bilinear-biquadratic SU{2) 
S=l model for d > 1. The phases are: AF: antiferromagnetic; SN: 
staggered-nematic, FM: ferromagnetic, and UN: uniform-nematic. 

In closing this section let us summarize the main steps to 
follow in order to obtain and classify the local order parame- 
ters. 

• Identify the group Qy\i = SU (D) associated to the hierar- 

chical language whose fundamental representation has 
the same dimension D as the local Hilbert space of the 
problem. The generators of this language exhaust all 
possible local order parameters. 

• Identify the group of global symmetries of the Hamiltonian 

Q which are direct products of local transformations. 

• Given that Q C Qy^i, then one can classify the generators 

of Qhi in the fundamental representation according to 
the irreps of Q. Each irrep leads to a different broken 
symmetry order parameter 

• Key: existence of a general set of SU {D) transformations. 
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IX. LATTICE GAUGE THEORIES AND THE QUANTUM 
LINK MODELS 

Wilson's lattice formulation of quantum field theory pro- 
vides a non-perturbative regularization to Euclidean func- 
tional integrals lf79il . The space-time lattice supplies the field 
theory with a cutoff such that loop integrations in perturbation 
theory yields finite results instead of divergences. 

In this way, a natural analogy emerges between quantum 
field theory on a d-dimensional lattice and classical statisti- 
cal mechanics in d + 1, which for a pure gauge theory (i.e., 
gauge theories without matter fields) amounts to consider the 
classical lattice action: 

S[u] = - Tr|ui^^Ui+/i^^u!^-^4,^| . (181) 

The classical complex parallel transporters ui transform as: 

= A-i(i)iii,^ A(i + A) (182) 

under local gauge transformations A(i) which leave the action 
S[u] invariant. The elementary parallel transporters are asso- 
ciated with an ordered pair of points in the lattice (links). For 
instance, the link (i, n) corresponds to the straight path from 
lattice site i to i + /t. 

Recently, Wiese and co-workers |8pll have further elabo- 
rated work started by Horn in 198 1 (81| on a new way of non- 
perturbative regularization of field theories which they named 
D-theory. In D-theory the Euclidean action in the standard 
Wilsonian formulation of a d-dimensional lattice field theory 
is replaced by a Hamilton operator H 

with quantum link operators Ui^^ on a d-dimensional lat- 
tice, which constitute generators of an algebra acting on a 
Hilbert space (e.g, a U{1) quantum link model can be real- 
ized with quantum links that satisfy an su{2) algebra). H 
commutes with the local generators of gauge transformations 
Gi = Dj A-i(j) C/i,^ niA(l) with A(i) = e-' ^O 
and the theory is defined through the quantum partition func- 
tion Z = Tr {e^*^^}, where the trace is taken on the Hilbert 
space. 

As emphasized in Ref. Isoll . D-theory is not a new set 
of field theories but another lattice regularization and quan- 
tization of the corresponding classical models. The main im- 
portant and attractive feature is the use of discrete quantized 
variables and finite Hilbert spaces. However, it turns out to 
be necessary to formulate the theory with an additional Eu- 
clidean dimension (i.e., in d + 1). In this regard dimensional 
reduction is a fundamental component to relate quantum link 
models to ordinary gauge theories. The dimensional reduction 
hypothesis relies on the existence of a massless phase ind+l, 
assumption that must be verified on a case-by-case basis. 

The existence of a quantum link model connecting lattice 
gauge theories to spin (or other algebraic) theories opens the 
possibility of formal connections between gauge theories of 



high-energy physics and strongly correlated problems of con- 
densed matter Indeed, since there is a connection between 
the "Gauge World" and the "Spin World" and we exhausted 
the connection between the "Spin World" and the "Particle 
World" the consequence is that one can find isomorphisms 
between the "Gauge World" and the "Particle World." In this 
way, for instance, one may look for the exact equivalent of 
confinement (in gauge theories) in the particle (condensed- 
matter) language. Even though some of these ideas were spec- 
ulated in the literature there was no formal rigorous relation 
established: The connection is our fundamental theorem. 

A. f/(l) Gauge Magnet 

For pedagogical purposes, it is convenient to illustrate the 
connection between gauge and condensed matter theories by 
considering the simplest U{1) gauge theory, usually called 
gauge magnet. Gauge magnets are gauge-invariant general- 
izations of the Heisenberg models. The corresponding classi- 
cal model has a J7(l) parallel transporter 

= e'*J- (184) 

associated to each link (j, fj) and the action S\'u\ is given by 
Eq. ( I18U . The dagger in this case denotes complex conjuga- 
tion. As it is shown by Chandrasekharan and Wiese 1,82,1. after 
the quantization process the uy^ fields become generators of 
a local su{2) algebra 

Ul, = Sr^ . (185) 

In other words, each link variable is replaced by a spin oper- 
ator Like in the case of quantum spin systems, this quanti- 
zation can be realized with any spin representation of SU{2). 
The simplest case corresponds to 5*= 1/2. In this case there 
are two possible states for each link variable which are de- 
noted by the two possible eigenvalues of S? ^ ~ ±1/2. The 
Hamiltonian of the quantum U{1) gauge magnet is obtained 
by replacing Eq. ( I185> into Eq. ( I183> 

and the generator of the U{1) gauge symmetry is 

G(i)=^(5f„^,^ + 5?^). (187) 

It is natural to ask what is the equivalent particle Hamilto- 
nian which is obtained by applying a spin-particle transforma- 
tion to Eq. ( I186t . We will only consider the simplest 5=1/2 
case for which Orland [82J obtained the exact solution in 2 + 1 
dimensions. The key observation made by Orland is that this 
spin model is equivalent to a gas of transversely oscillating 
fermionic strings. The string tension can be obtained from the 
exact solution and it was shown that the fermionic charges are 
confined. To obtain this result we need to introduce a dual 
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transformation which replaces the links by sites: (i, ly) r, 
(i, fi) r + X, and {i + 0, fi) r + y. This transformation 
is illustrated in Fig.^J where we can see that the dual lattice 
is divided into a checkerboard pattern. This means that a pla- 
quette with vertices r = {ri,r2), r + x, r + x + y, and r + ij 
is black if ri + r2 is even and white if ri + r2 is odd. The 
sublattice of points r such that ri + r2 is even will be denoted 
by A. 




FIG. 12: Transformation from the direct square lattice, with lattice 
sites i and links (i, (i, v), to the dual checkerboard lattice, with 
sites r and links (r, x) and (r, y). 

The expression for Hgrn after the dual transformation is 

reA 

+ ^r+x+y^r+x^r ^r+y ' (188) 

It is clear from this expression that the black plaquettes are 
active and the white ones are passive. Equation ([^^ implies 
that the total charge in each white (passive) plaquette is con- 
served. In absence of background sources, the physical states 
are those which are annihilated by the 'Gauss law' operator 
G(i). In the particle language, this condition means that there 
are two particles on each white plaquette (see for instance the 
state illustrated in Fig.ll3>. 

We have seen in section lVT!! that a spin 5=1/2 can be trans- 
formed into a spinless hard-core anyon through the general- 
ized JW mapping 

S- = (At)^ar, 

= a^flr - \ ■ (189) 

For 9 = 0, the anyons become hard-core bosons and for 6 = n 
they are spinless fermions. By replacing these expressions 



into Eq. ( I188t we get the following Hamiltonian representing 
interacting hard-core anyons 

reA 

+ a^+x^v+x+y^l+yal ■ (190) 

Notice that the transmutators and (A'^)^ do not appear 
in the Hamiltonian. This means that Hgrn is invariant under 
transmutation of the statistics. To physically understand the 
origin of this invariance we just need to realize that the dy- 
namics imposed by H^rn only allows the motion of particles 
along the horizontal axis within a passive plaquette (see Fig. 
I13> . Therefore, the original ordering of the particles is pre- 
served and the statistics turns to be irrelevant. 

The exact ground states of this model l83ll correspond to 
parallel strings which are aligned in the y direction (see Fig. 
I13> . The model can be exactly solved because the quantum 
fluctuations of each string are described by an effective one- 
dimensional 5=1/2 XY model l83ll . Notice that this solu- 
tion relates to the new paradigm in strongly correlated mat- 
ter where the appearance of local inhomogeneous (stripe-like) 
structures seems to be a common feature of many different 
physical systems. 

There are other possible connections between the gauge 
magnet and condensed matter theories which correspond to 
different changes of language. If each link is associated to 
a singlet bond state (Cooper pair), Orland l84ll has shown 
that the 5*= 1/2 gauge magnet is mapped onto the Rokhsar- 
Kivelson model L85.1 (without the diagonal term) if the new 
'Gauss Law': [G{i) + 2]\^) = is imposed. In this way 
we see how the formal connections between condensed mat- 
ter and lattice gauge theories can be exploited to predict new 
physical behaviors in one or the other field. 



String phase 



FIG. 13: Schematic picture of a ground state of the S=l/2 gauge 
magnet in the particle language. Quantum fluctuations (QF) corre- 
spond to the horizontal motion of vertical pairs in the active (black) 
plaquettes. 



X. QUANTUM INFORMATION AND COMPUTATION 

A new challenge in information theory and computer sci- 
ence has recently emerged as the result of applying the fun- 
damental laws of quantum mechanics and using the quantum 



39 



LA-UR-OO-XXXX 



February 6, 2008 



Advances in Physics 



effects to advantage. This new set of ideas comprise what 
is known as "Theory of Quantum Computation and Quantum 
Information" and has as a major objective to process infor- 
mation in a way that exceeds the capabilities of classical in- 
formation LlSiil . The device that performs the manipulation of 
information is named quantum computer and the standard unit 
of information is the qubit (i.e., a two-level system). The close 
relationship between information processing and the physical 
phenomena leading to it is perhaps the most remarkable aspect 
of this new paradigm. Since information can be represented 
in many different physical forms, and easily converted from 
one form to another without changing its meaning, quantum 
information represents a new abstract archetype for informa- 
tion processing independent of the precise implementation of 
the quantum computer, only requiring at least one physical 
representation to be useful. 

A key fundamental concept in information theory is the re- 
alization (T3l that a model of computation is intimately con- 
nected to a physical system through a closed operator alge- 
bra. In other words, each physical system is associated to a 
certain language (e.g., spin S=l/2) and thus to an algebra re- 
alizing it (e.g., Pauli algebra), and that particular algebra may 
become a possible model of computation. An immediate con- 
sequence is that an arbitrary physical system can be simulated 
by another physical system (e.g., a quantum computer) when- 
ever there exists an isomorphic mapping between the different 
operator algebras representing the systems tl61. A very sim- 
ple example is provided in Refs. jj^ and l86ll . where it is 
shown how to simulate aid impurity Anderson model using a 
quantum computer based on an array of spins S=l/2. Another 
resource of quantum information that can help us understand 
complex phenomena in strongly correlated matter is the no- 
tion and measures of entanglement. 

It is very clear the power that our fundamental theorem has 
by providing the formal connections (isomorphisms) between 
the different languages of nature. Therefore, the implications 
for quantum information and computation are rather obvi- 
ous, namely that one can identify quantum resources and de- 
fine convenient models of computation, or imitate an arbitrary 
quantum phenomena with a given quantum computer given 
the appropriate dictionaries to translate nature's language to 
the machine language. In this way, one can recognize three 
complementary areas where the concepts of language and dic- 
tionaries developed in the present manuscript are of particular 
relevance, areas we wiU expand on in the following. 

Models of computation and quantum resources 

A model of computation consists of an algebra of opera- 
tors, a set of controllable Hamiltonians and measurable 
observables, and an initial state of the physical system. 
The set of controllable Hamiltonians must be universal 
in the sense that any unitary operation can be performed 
with such a set. The standard model of computation is a 
particular example which uses the Pauli operator alge- 
bra and the qubit as basic unit. Another possible model 
is the fermion model fl6ll which is isomorphically re- 
lated to the standard one through the JW mapping and, 
thus, is equivalent. The main point is that the choice of 



model depends upon the nature of the available physical 
resources and their quantum control. For example, in a 
liquid NMR quantum computer the nuclear spins of the 
molecules (5*= 1/2) are the units which can be controlled 
and the standard model is the appropriate one. 

The control of quantum mechanical systems is ham- 
pered by quantum noise and decoherence and, there- 
fore, identifying static quantum resources for informa- 
tion processing is a great challenge. Given a physical 
information processing device there will be a language 
which will the most natural one for the elementary static 
resources of the device. For instance, if we could con- 
trol "'He atoms at the quantum level it would be natural 
to consider a hard-core boson model of computation. 
The importance of our fundamental theorem is that once 
one identifies the best quantum resources it allows us to 
build the model of computation accordingly. 

Simulation of physical phenomena 

Physical phenomena can be simulated or imitated by 
a quantum network |87] with the help of a quantum 
computer Imitation is realized through a quantum al- 
gorithm which consists of a quantum network with a 
means to repeat blocks of instructions. A quantum net- 
work is defined by a sequence of universal gates (uni- 
tary operations), applied to the system for the purpose 
of information processing, and measurements in a fixed 
temporal order The measurement operation is mostly 
needed to classically access information about the state 
of the system. Every matrix which represents a re- 
versible operation on quantum states can be expressed 
as a product of the one and two-qubit gates, and the 
minimum set needed to represent any such matrices is 
called a universal set of gates. 

When trying to simulate a problem using quantum in- 
formation processing, an important issue is to determine 
how many physical resources are needed for the solu- 
tion. The main resources are quantum space, the num- 
ber of qubits needed, and quantum time, the number of 
quantum gates required. The accounting of algorithmic 
resources forms the foundations of quantum complexity 
theory. One of the objectives in quantum information 
theory is to accomplish imitation efficiently, i.e, with 
polynomial complexity, and the hope is that quantum 
imitation is more efficient (i.e., needs less resources) 
than classical imitation. There are examples that sup- 
port such hope (e.g., fermion simulations with poly- 
nomially bounded statistical errors ITS IstIi ). although 
there is no general proof that indicates the superiority of 
quantum over classical imitations, regarding efficiency. 
Indeed, there is, so far, no efficient quantum algorithm 
that can determine the ground state (or, in general, the 
spectrum) of a given Hermitian operator ISTII . despite 
occasional claims. It is known that the ability to re- 
solve this question leads to efficient algorithms for NP- 
complete problems like the traveUng salesman conun- 
drum. 
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A very important observation, in connection with the 
notion of efficiency, is a corollary of our fundamental 
theorem: Given two languages, the generators of one 
of them can be written as a polynomial function, with 
polynomial complexity in the number of modes or re- 
sources, of the generators of the other and vice versa. 
This result implies that the important algorithmic step 
of translation from the language of the system to be 
imitated to the machine language does not change the 
complexity of the quantum space and time. 

Certainly, a general purpose quantum computer is not 
the only device that allows simulation of physical phe- 
nomena in nature (with its many languages). Imitation 
can also be achieved in a conceptually different manner 
using a quantum simulator The main distinction is the 
lack of universality of the latter An example of a quan- 
tum simulator is an optical lattice |88] which is specif- 
ically designed to imitate a given physical Hamiltonian 
and where there is limited quantum control. The possi- 
bility of control and tunability of the interactions of the 
elementary constituents offers the potential to design 
new states of matter l60Ll88il1 . This is of particular rel- 
evance in strongly correlated matter where these quan- 
tum simulators furnish the benchmark to test theories 
and approximations L6Q1I . Again, the importance of the 
languages and dictionaries developed in this manuscript 
is clear and concrete. 

Quantum information measures 

Entanglement, a word that Schrodinger coined to dis- 
tinguish quantum from classical mechanics, is that 
bizarre feature of composite quantum systems that led 
to so much controversy in the past (like quantum non- 
locality in the Einstein-Podolsky-Rosen gedanken ex- 
periment). Essentially, entanglement is a quantum 
property whereby a pure state of a composite quantum 
system may cease to be determined by the states of its 
constituent subsystems. Entangled pure states are those 
that have mixed subsystem states. An independent con- 
cept is the notion of separability which refers to the 
property of a quantum state of a composite system of 
being able to be written as a direct product of states 
of its component subsystems. The theory of entangle- 
ment, and its generalizations, is currently under devel- 
opment fTM l89il . In Ref . l89il entanglement is viewed 
as an observer-dependent concept, whose properties are 
determined by the expectations of a distinguished sub- 
space of observables of the system of interest (exper- 
imental access), without reference to a preferred sub- 
system decomposition. The standard notion of entan- 
glement is recovered when these means are limited to 
local observables acting on subsystems. A tremendous 
effort is put in trying to understand the properties of en- 
tanglement that can be used as a resource in quantum 
information. 

In a sense, that must be defined more precisely Is^ . 
the notion of entanglement is a relative of the notion 
of quantum correlations and a very relevant question 



is whether one can construct useful measures of entan- 
glement to understand the emergence of complex phe- 
nomena in strongly correlated matter A quantum phase 
transition involves a qualitative change in the correla- 
tions of the ground state of the system as a result of 
tuning parameters of its Hamiltonian. In some cases an 
order parameter is associated to the transition, in others 
a topological order It is intuitively expected that this 
change is also associated to a change in the nature of 
entanglement. Therefore, quantifying and classifying 
entanglement is very important to characterize a quan- 
tum phase transition: Can measures of entanglement 
distinguish between broken and non-broken 1 2] symme- 
try phase transitions? 



XI. SUMMARY 

We have introduced an algebraic framework for interacting 
quantum systems to study complex phenomena characterized 
by the coexistence and competition of various ordered states. 
We argued that symmetry, and topology, are key guiding prin- 
ciples behind such complex emergent behavior Emphasis has 
been made in developing a systematic mathematical structure 
that allows one to attack these problems within a single unify- 
ing approach. 

The core result of the paper, from which all other results 
follow, is the proof of a fundamental theorem that permits 
to connect the various operator languages used in the de- 
scription of the properties of physical systems. This theo- 
rem together with the notion of transmutation of statistics pro- 
vide the tools necessary to unifying the quantum description 
of matter To formalize this unification we needed to rig- 
orously define the concepts of language and dictionary: To 
model a particular physical phenomena we commonly iden- 
tify the main degrees of freedom of the problem and associate 
to them certain operators. One can furnish to the resulting set 
of operators (that we call language) with an algebraic struc- 
ture and ask whether two different languages have something 
in common. The fundamental theorem tells us that two lan- 
guages can be connected whenever the dimension of their lo- 
cal Hilbert spaces are equal. We expanded the notion of local 
Hilbert space to embrace different Hilbert space decomposi- 
tions (we saw, for instance, how to map the Hilbert space of 
a bond to a site). The resulting one-to-one language map- 
pings we named dictionaries (a traditional example of which 
is the Jordan- Wigner mapping). In the course of the presen- 
tation we showed, through example, many different dictionar- 
ies relating diverse operator languages. In this way we de- 
fined universality of behavior as an equivalence relation be- 
tween seemingly different physical phenomena which share 
exactly the same underlying mathematical structure as a result 
of one-to-one language mappings (for example, the spin ne- 
matic order and Bose-Einstein condensation of flavored hard- 
core bosons). Out of the many languages one can use to de- 
scribe a given physical problem there is a class, we named 
hierarchical language, which has the advantage that any local 
operator can be expressed as a linear combination of its gen- 
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erators. In this way, hierarchical languages provide the tools 
necessary to classify order parameters. 

A new, formally developed, notion is the idea of emer- 
gent symmetry, i.e., the fact that new symmetries not realized 
in the Hamiltonian describing the system can emerge at low 
energies. From a broader perspective, an emergent symme- 
try is a non-local transformation which commutes with the 
system Hamiltonian and becomes local when restricted to a 
low-energy subspace of the Hilbert space. In some instances 
like the quasi-exactlv solvable problems, to which the Id t- 
Jz model belongs 01, or the family of spin Hamiltonians for 
which the ground state is a product of spin singlets (e.g., the 
Majumdar-Ghosh model), the low-energy effective theory can 
be derived in an exact way and the emergent symmetry is ex- 
act. In other cases, the emergent symmetry is approximate, 
however, it still provides a guiding principle to identify the 
nature of the ground state and its low-energy excitations. In 
the same way that rigidity is related to a broken symmetry 
phenomenon, localization (in real or any other space) is con- 
nected to the emergence of a gauge symmetry. Indeed, the 
metal-insulator transition may be interpreted as a discontinous 
change in the group of emergent symmetries characterizing 
the two phases. 

Figure ^5 summarizes the spirit and fundamental concepts 
that emerge from our fundamental theorem. 



There are several reasons for our algebraic framework to 
constitute a powerful method to study complex phenomena 
in interacting quantum systems. Most importantly: To con- 
nect seemingly unrelated physical phenomena (e.g., high-Tc 
or heavy fermions and quantum spin theories); identify the 
general symmetry principles behind complex phase diagrams; 
unveil hidden symmetries (and associated order parameters) 
to explore new states of matter; obtain exact solutions of rele- 
vant physical models that display complex ordering at certain 
points in Hamiltonian space; and find new approximations 
which do not favor any of the competing interactions. The 
power of the present approach is reflected in the unlimited 
number of potential applications where they could be used, 
and ranges from condensed matter and statistical mechanics 
to lattice gauge theories and quantum information and com- 
putation. We hope and envision this formalism will find fu- 
ture applications in contexts different than the ones we have 
developed here. 
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